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Abstract: We study the five-body decays fi~ e~e^e~u^Ue and t~ -a 
for £' = e, within the Standard Model (SM) and in a general effective field theory 
description of the weak interactions at low energies. We compnte the branching ratios 
and compare onr resnlts with two previons, mntnally discrepant, SM calcnlations. By as- 
snming a general strnctnre for the weak cnrrents we derive the expressions for the energy 
and angnlar distribntions of the three charged leptons when the decaying lepton is polar¬ 
ized, which will be nsefnl in precise tests of the weak charged cnrrent at Belle II. In these 
decays, leptonic T-odd correlations in triple prodncts of spin and momenta -which may 
signal time reversal violation in the leptonic sector- are snppressed by the tiny nentrino 
masses. Therefore, a measnrement of snch T-violating observables wonld be associated to 
nentrinoless lepton flavor violating (LFV) decays, where this effect is not extremely snp¬ 
pressed. We also stndy the backgronnds that the SM five-lepton lepton decays constitnte 
to searches of LFV L~ -A decays. Searches at high valnes of the invariant mass 

of the pair look the most convenient way to overcome the backgronnd. 
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1 Introduction 

Five-body decays of leptons, L~ —)■ £~f~^£'~ULi^e with L = t, fi and £, £' = e, /r, are 
allowed processes in the SM, suppressed by a factor with respect to the lowest order 
three-body L~ —)■ £~viV(^ decays. Their study is interesting because they can provide an 
important source of background for searches of lepton flavor violating /i —)■ ee+e“ [1] and 
r ££'^£'~ decays due to undetected neutrinos. Analyses of these decays are also aimed 
at searches for sterile neutrinos and dark photons [2, 3]. 

Besides, these processes allow for stringent tests of the weak charged current and, 
in particular of its Lorentz structure and possible violations of lepton universality (for 
a recent review on related studies in tan lepton decays see [4]). In the closely related 
radiative L~ —)■ £“z/£z7^7 decays, the photon carries information about the outgoing lepton 
polarization and, as a result, two additional Michel-like parameters [5, 6] can be extracted 
with precision, in contrast to the non-radiative channel. For the tan lepton case, they 
shall be measured with comparable precision to the corresponding muon decays [7] using 
Belle-II data [8] thanks to the expected Belle-II statistics and performance [9]. These 
bright experimental prospects demand a corresponding effort on the theory and Monte 
Carlo [10] side. In this respect, let us mention the recent papers [11, 12] studying radiative 
/i and r leptonic decays at next-to-leading order and polarized r —)■ 3£ decays (focusing 
on angular correlations of new physics operators), respectively. 

In view of these forthcoming good quality data sets, it is timely to attempt to improve 
the current description of these decays. Among them, we shall include for the first time 
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Figure 1. Feynman diagrams for five-body decays of tan leptons. For identical leptons = i) 
in the final state, two additional diagrams corresponding to the exchange pi -f-)- p 2 should be 
considered. 

mass effects of the daughter charged leptons in decays of polarized particles. We will 
also address the accuracy reached, re-evaluating the branching ratios of different decay 
channels and comparing with existing (and conflicting) predictions. We hrst consider 
the contribution of lepton-flavor conserving charged weak interactions to these processes 
within an effective held theory framework, so as to test with precision the SM V — A 
universal structure of the W current through Michel-type parameters both in the spin- 
independent and spin-dependent cases. 

T-odd correlations in these leptonic decays seem to be absent in the SM and in its 
most general, lepton-havor conserving, effective held theory extension here considered. 
Very light neutrino masses (and thus hxed helicity, to an excellent degree of accuracy) 
may lead to negligible T-odd spin-momenta correlations in these decays however, no 
gain is obtained by relaxing the lepton-havor conserving requirement as it is shown by 
considering a simple example of a new physics extension with Z' gauge-boson. Therefore, 
an eventual observation of these T-violating correlations may provide an indirect signal 
of a non-conventional source of CP violation. This one, in turn, would contribute to 
generate a net baryon asymmetry in early stages of the universe [13] through the anomaly- 
free character of the B — L accidental symmetry of the SM [14, 15], with fundamental 
implications for the enormous matter-antimatter asymmetry of the universe and cosmic 
evolution. 

Previous studies of the branching ratios for hve-body decays of unpolarized leptons 
within the SM have been reported in [16-23]. While diherent calculations of fi decays 
agree within numerical uncertainties, the results of two theoretical calculations which 
consider tau decays [20, 21] differ signihcantly. Numerical studies of the angular and 
energy distributions of charged leptons in p decays have been reported in [18] by assuming 
the most general form of lepton-havor conserving weak currents including (axial-)vector, 
(pseudo)scalar and tensor couplings. In the present study, we compute the branching 

^ T-odd correlations induced by genuine T violation would provide indirect signals of CP violation in 
the leptonic sector. 
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ratios of five-body leptonic decays of tau and muon leptons using an integration method for 
the phase-space that differs from previous studies. We also provide analytical expressions 
for the angular and energy distribution of charged leptons in the case of the SM couplings 
for decays of polarized leptons keeping the whole charged lepton mass dependence but 
always neglecting neutrino masses. We also compute these distributions in the case of the 
most general Lorentz structure for the charged weak currents and provide their analytic 
expressions, in terms of the Michel parameters introduced in Ref. [17], to describe five- 
body decays in the polarized case keeping hnite daughter lepton masses for the hrst time. 
Our expressions are aimed to allow stringent tests of the SM charged weak current by 
the experimental collaborations. Then, we provide an example of a new physics (LFV) 
extension which shows how T-odd spin-momenta correlations are absent, as in the case of 
the charged weak current with the most general Lorentz structure, in the limit of massless 
neutrinos. The observation of T-violation in these observables would therefore hint at the 
corresponding neutrinoless LFV processes, which strengths the case for these searches 
in leptonic muon and tau decays. Finally, we consider in detail the backgrounds that 
five-lepton lepton decays constitute for searches of the corresponding neutrinoless LFV 
processes and summarize our conclusions. The explicit expressions of all contributions to 
the observables that we have computed are collected in the Appendices. 

2 Five-lepton channels within the SM 

In the SM, the hve-body decays of muons and tans can be produced, at lowest order, from 
a lepton-pair conversion of a virtual photon emitted from the usual leptonic decay as 
shown in Figure 1. For identical charged leptons in the hnal state, two other diagrams 
contribute such that the decay amplitude becomes anti-symmetric under their exchange. 

2.1 Notation and kinematics 

The convention of momenta for the illustrative case of r lepton decays, is 

r"(Q) r{pi) i'~{p 2 ) MPi) ^riPb) ■ ( 2 . 1 ) 

Similarly, the masses of the decaying and hnal-state charged leptons will be denoted by 
M, mi, m (m 2 = m 3 = m), respectively. In the SM, the charged weak currents have a 
V — A structure and their strength is encoded in the Fermi constant Gp- 
For non-identical fermions (F 7 ^ £), the decay amplitude becomes: 

M-sm = /- ^ (Ml -I- M. 2 Y , (2.2) 

^Lepton-pair conversion from Z and Higgs bosons replacing the photon are also possible, but are 
extremely suppressed with respect to electromagnetic conversion, as shown in Figure 1 for the r lepton 
case. 
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with 


-Ml = M(P5)7a(l - 75 ) ^_m ) “ 75)^^(P4) , (2.3) 

-M^ = m(p5)7“(1 - 75)m(Q) • h(pi)7'' { , ,] -1 7a(l - 75)w(P4), (2.4) 

\pi+ i-mij 

where = u{p 2 )l^v{p^) is the electromagnetic vertex current and (? = P 2 +P 3 is the mo¬ 
mentum of the virtual photon {q^L^ = 0 due to conservation of electromagnetic current). 
For identical leptons (£' = £), the decay amplitude becomes M. = M.i + — + M. 4 )i 

where the last two terms are obtained from the hrst two by the exchange pi ^ P 2 - 

Owing to the masslessness of neutrinos, the generic form of the unpolarized squared 
matrix element can be written as: 

= = ( 2 . 8 ) 

pols 

This tensor structure is preserved in the case of polarized decaying particle. 


2.2 Branching ratio 

The decay width in the rest frame Q = (M, 0) of the decaying particle is ^ 

^-w/n^wF8‘(Q-|:p. 


( 2 . 6 ) 


The squared (unpolarized) amplitude |A4P of the hve-body decay depends upon 8 inde¬ 
pendent kinematical variables. According to Ref. [24], we choose them to be: 

Sl = (Q-Pl)^ S2 = (Q - Pi - P2)^ S 3 = (Q-pi-p2-P3)^ 

UI = {Q-P2f, U2 = {Q-P3f, U3 = {Q-Pif, 

t2 = {Q -P2-P3f , h = {Q -P2-P3-PAf ■ (2.7) 


For the reader’s convenience, we quote the expression of the relevant scalar products in 
terms of our set of kinematical invariants given by eq. (2.7). 


Q-Pi = 


-|- m\ — Si 


T 


, Q-P2 = 


— ui 


T 


, Q-P3 = 


_ M^ -I- — U 2 


T 


P „ _ — U3 „ _ S2 — Si — Ml -|- _ S3 — S2 — ^2 + Wi 

V P4 — - 2 - ’ Pi P 2 — - — 2^ - ’ Pi P3 — -^ 


to — — to 

Pi‘P a = ^^^ 


U — u^ — U2 + 

P2-P3 = — - 2 - 


— M 3 — ^2 + ^3 — 2p3 • P4 

P2 ■ Pa = -^^- 


( 2 . 8 ) 


factor of 1/2 in the right hand-side needs to be added when dealing with identical particles in the 
final state. 
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The very long expression for ■ p/^ is qnoted in Appendix C, for completeness. We have 
verified that our expression for ps ■ p^ agrees with eqs. (A.6) to (A.8) in Ref. [25] with 
suitable replacements for the squared masses {i = 1,.., 5). 

Our method for integrating the phase-space to get the branching fraction is different 
from the one used in Ref. [20]. We integrate the phase-space directly over the eight 
independent kinematical variables mentioned above. In contrast, in Ref. [20] the phase- 
space is first (partially) integrated over the momenta of the two neutrinos using the 
covariance properties of the tensor integral {P = Q — pi — P 2 — Ps), namely 


48(27r)5 


I (P-P.- K) 

^ P^ 9ap + B PaPfS , 


(2.9) 


with B = 2A = 1/ (24(27r)®). Then, the integration over the remaining variables is 
carried out. Therefore, the phase-space integration using our special choice of kinematics 
will allow to obtain an independent verification of previous calculations in [20]. 

The matrix element squared and summed over all lepton polarizations reads 


|W1|2 = e^Gi 


Tii+ T22 + T1221 + T33 + ^44 — ^1331 ~ ^1441 — ^2332 “ ^2442 + ^3443 


( 2 . 10 ) 


Tii corresponds to the contribution of MiAij, while Tijji to that of MiAij -j- AijAij. 
The subindexes 1 and 2 stand for Figs. 1 (a) and (b), respectively, while diags. 3 and 
4 are obtained -in this order- by exchanging identical fermions in diagrams la and lb. 
Subindexes in eq. (2.10) stand for the corresponding contributions of different diagrams. 
Analogous notation will be used throughout. For convenience, we define the reduced 
amplitudes (ta, tijji) taking out common factors 


Til = 

Ti221 = 


321 


11 


[T>i(P2)T>2(P2)]^ 
32 ti221 

[D,ip2)]W2{p2)D 

32^44 


T44 — 


Ti441 — 


[Ti(pi)T]2 
32 ti44i 


T>i(Pi)T>i(P2)T>2(P2)T> 
32 ^2442 


T2442 — 


D,{p,)D,{p2)D^ 


T 22 = 
T33 = 
Ti331 = 
T2332 = 
T3443 = 


321 


22 


[Di{P2 )D? ’ 
32 ^33 


( 2 . 11 ) 


[D,{p,)D2{p,W ’ 

32 ti33i 


T>i(Pi)T>i(P2)T>2(Pi)T)2(P2) ’ 
32 ^2332 

T>i(Pi)T>i(P2)T>2(pi)T) ’ 

32 ^3443 

[D,{p^)]W2{p,)D • 


Explicit expressions for the tu and tijji are given in Appendix A. 
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In eq. (2.11) we have introduced additional variables to make expressions shorter. 
These are {i = 1,2) 


Di{pi) = + Pi ■ , ( 2 . 12 ) 

D2{pi) = m'^ +pi-p3-Q -pi-Q -ps, 

D = + Pi ■P2 +Pi ■pz+P2-Pi.- 

Note that these definitions are true for both decay modes, without identical particles and 
with identical particles (mi = m). 


Ghannel 

Ref. [20] 

Ref. [21] 

This work 

PDG [27] 

BR(r“ —)■ e“e’''e“t'ei^r) X10^ 
BR(r“ e“/i’''/i“i^eG)xl0^ 
BR(r“ —)■ p,~e~^e~iifj,i'r)^^0^ 
BR(r“ — 

BR(yU“ e“e+e“t'et'At)xl0^ 

4.15 ±0.06 

1.257 ±0.003 

1.97 ±0.02 

1.190 ±0.002 

3.60 ±0.02 

4.457 ±0.006 

1.347 ±0.002 

2.089 ±0.003 

1.276 ±0.004 

3.605 ±0.005 

4.21 ±0.01 

1.247 ±0.001 

1.984 ± 0.004 

1.183 ±0.001 

3.597 ±0.002 

2.8 ± 1.5 

< 3.6 

3.4 ±0.4 


Table 1. Branching ratios for the five-body decays of r and p leptons with quoted error bar 
obtained from accuracy of numerical integration. Some of the previous calculations are shown, 
for comparison, in the second and third columns. An additional ±0.17% uncertainty owing to the 
current precision of the r lifetime measurement [27] should be added to our result. Experimental 
data are scarce, with large error bars, but still consistent with the SM predictions. 


Our results for the branching ratios are shown in Table 1 and are compared to previous 
results and available experimental data. The quoted uncertainties arise from integration as 
given by VEGAS [26], the numerical code used in our numerical calculations (as well as in 
Ref. [20]). Table 1 shows that our results are in good agreement with those of Ref. [20], and 
substantially differ in some cases from those of Ref. [21]. Only tiny differences might be 
attributed to the use of different inputs. Unfortunately, we cannot read from those papers 
the exact values of the employed inputs. We use the values of masses -and particularly 
the tan mass, (1776.82 ± 0.16) MeV- quoted in the 2014 edition of the PDG [27], while 
Mr of Ref. [20] (and [21]) presumably corresponds to the GLEO 1993 measurement [28] 
(1777.8±1.8) MeV. Analogously, the current tan lifetime, (2.903±0.005)-10“^^ s, although 
consistent with the ALEPH 1992 measurement [29], (2.91 ± 0.14) ■ 10“^^ s, is by far more 
precise. Our updated values may explain the small differences between our results and 
those of Ref. [20], but not with those of Ref. [21]. In particular, focusing in the comparison 
of our numerical results with those of Ref. [20], we notice larger differences in the modes 
with two or three muons in the final state. Taking into account that different tan lifetime 
input results in an overall shift for all modes, and that modes with heavier daughter 
leptons are more sensitive to the tan mass owing to kinematics, we can attribute the 
small differences to the used tan mass inputs. 


6 









In Fig. 2 we plot the normalized differential decay width versus Ei = E^- for all five- 
lepton r decay channels. As expected, no dynamical structures can be seen and only the 
kinematical enhancement in the 3e channel for low electron energy is noticeable. In Fig. 3 
we plot the normalized differential decay width versus for the considered 

decays. In this case, all channels tend to peak at low values of because of the low 

virtuality of the photon near threshold. 


t 

> 

CD 

O 




Ej- (GeV) 


Figure 2. Differential decay width (normalized to the partial decay width) versus E^- (in modes 
without identical particles it corresponds to the energy of the different charged lepton) for all 
five-lepton r decay channels. 


Our precise predictions can also be useful to provide an independent test of an 
anomaly recently reported in tan leptonic decays. In particular, BaBar’s measurement 
of r —)■ [30] differs by 3.5cr from the SM prediction [11], given at NLO. On the 

contrary, the agreement on the muon mode is at the 1 a level [30]. An independent study 
of these decays with Belle-I data could help to settle this issue. Given our precise predic¬ 
tion for the r SeUei^r decays and its branching ratio (~ 4 ■ 10“®), measurable with hrst 
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generation B-factories BaBar and Belle-I data, another interesting check of the anomaly 
would come from analysing the former process. 

2.3 Spin-momenta correlation in decays of polarized leptons 

The polarization of the decaying particle is introduced by replacing 

+ M) —t -(^ + M)(l + 75 /) 

appropriately in the unpolarized squared amplitude. The polarization four-vector 
satisfies the properties Q • s = 0 and = — 1. In the rest frame of the decaying particle, 
s = (0,s). 

Considering the polarization of the decaying particle does not change the dependence 
of the squared amplitude upon the neutrino four-momenta ^ 4 , 5 , as compared to the unpo¬ 
larized case, eq. (2.5). Therefore, the partial integration of the phase-space over neutrino 


8 











four-momenta is given by Eq. (2.9), and the differential decay rate by: 


dT, 

#pid3p2(iSpj 3. 2'V>MEtE2E3 ’ 


(2.13) 


where we have added the upper-index s on 7)^^^ to make explicit the dependence on the 
tau polarization. 

Thus, after integration over neutrino four-momenta the differential rate depends 
only upon the product of charged particle four-momenta {Q,Pi,P 2 ,P 3 ) and the polarization 
four-vector s. In the rest frame of the decaying particle, the numerator in the right hand 
side (r.h.s.) of Eq. (2.13) can be written as 




L Pi ■ s — Gi p2 • s — G2 Ps ■ s 


(2.14) 


In the above equation, the coefficients F, L, Gi and G 2 depend only upon dot products 
of charged leptons momenta. 

In order to cast our results in a more compact way we will rewrite the previous result 
using the most commonly used phase-space integration variables 


dr, 


dxidVtidx2dVt2dx‘j,dVt^ 


lll^|Pl||P2||P3| 

3.22ijrio 




(2,16) 


which also bear a closer relation to the measured observables. In the previous equation, 
the reduced dimensionless variables Xi = 2Ei/M {i = 1,2,3) were introduced. 

In this way, the matrix element squared, summed over hnal-state lepton polarizations 
and integrated over the neutrino phase space reads 


TUr^{P) = e-Gj. 


T„ + T, 


22 


-1221 


+ ^33 + Ti 


44 


-1331 


-1441 


(2332 


= e^GlF 


F-'-SM , 


T2442 -l- T3443 
(2.16) 


where we made explicit the relative signs due to exchange of identical fermions For 
later convenience we will use the reduced and tijji amplitudes which are dehned as ta 
and Ujji in eq. (2.11). 

^As requested by the Belle Collaboration, we have also provided them with expressions for the cor¬ 
responding spin (in)dependent form factors without integrating over neutrino four-momenta. This is 
intended for direct implementation in the TAUOLA Monte Carlo Generator [31-34] for use of the col¬ 
laboration. We are not including these expressions in this paper, but will send them upon eventual 
request. 

®This notation corresponds to Belle’s conventions and is useful to display relationships among Gi 
and G 2 in the case without identical particles. Corresponding symmetries are found between L and G 2 
when there are indistinguishable fermions. F can be checked with earlier computations, i.e. [20], for 
the unpolarized case with the appropriate changes of kinematical invariants. We present, however, our 
complete expressions to avoid the reader these uncomfortable translations. 

®In case all particles are distinguishable, only the first three terms contribute in eq. (2.16). 
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Taking advantage of the decomposition in eqs. (2.14) and (2.16) we can split the 
different contributions to the spin-(in)dependent form factors as 

F = Fii + F22 + F122I + -^33 + -^44 — -^1331 ~ -^1441 — -^2332 “ -^2442 + -^3443 ; (2-17) 

and analogously for Gi, G 2 and L. We can further introduce the reduced form factors 
/j Qii 92 and I functions by factoring out the same common coefficients as in eq. ( 2 . 11 ). 
The explicit expressions of these reduced form factors are given in Appendix B. Trace 
identities [35, 36] have been used to simplify the calculations. 


3 Effective field theory analysis 


In the previous section we considered in detail the SM predictions for five-body leptonic 
decay modes of polarized and unpolarized /i and r leptons. Here, we will generalize this 
approach within an effective quantum field theory description of the weak charged current 
at low energies, i.e., much smaller than the electroweak scale. We focus on decays of 
polarized /r, r leptons as the different spin and charged lepton momenta correlations may 
be useful to study the effects of New Physics described by new operators and couplings 
of the effective weak Hamiltonian. 

The most general local, derivative-free, lepton-number conserving Lagrangian describ¬ 
ing four-lepton interactions consistent with locality and Lorentz symmetry can be written 
[5, 6 , 17, 18, 37-39] 


C 


AGii' 

W 


i,X,p 




(3.1) 


with i = S', y, T; T"^ = /, = 7 ^, T^ = /\/2 labelling the Lorentz structure of weak 

currents and X, p = L, R the chiralities of the charged leptons (.^ and k, the chiralities of 
neutrinos, are hxed by A and p once a particular Lorentz structure T* is chosen). There 
are 10 independent complex coefficients (4 scalar, 4 vector plus the two tensors-type 
couplings g"[j^ and which altogether give rise to 19 independent real couplings once 
an unphysical global phase is removed The global factor Gni is hxed by the width of 
t~ decay, in such a way that [17] 

1 = I X] IS'fpP + IS'ApP + ^{\9 rl \'^ + \9 lr \'^) ■ (3-2) 

A,p A,p 

The decay amplitude for hve-lepton channels in r and p decays in the general case will 
receive contributions from the different terms of the Lagrangian (3.1). Those contributions 

^The g\p coefficients, which parametrize beyond the SM effects at low energies in the weak charged 
current, can be related to the Wilson coefficients of the effective Lagrangian at the electroweak scale 
[40] (see [41] for an updated discussion in the context of tau lepton decays). This allows to complement 
low-energy results with high-energy LEP/LHC searches. 
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to will be denoted as follows 


T:^,I'^i>{P) = e*|G„f [is 


S \2rriS I I^V 
-t ; 


LL 

T* 


LL + \9lL 
ST I S 


I \„S 

I -L I 


LL + \9rl 
y* rpsv 


2t^s I i^y 
I -L i 


RL + \9rl 
S „V*rrSV 


IrpV 
^ RL 


+ \9_ 


T 

RL 


2rpT 

^RL 


-LSPo I I 5 V * rpHV i S V * rpHV i V T * rjiVT 1 

s--^^\9rl9rl ^rlrl “t 9ll9rr -‘-llrr “t 9lr9rl -‘-lrrl “t 9lr9rl -‘-lrrl J 


-\-L R 


(3.3) 


in such a way that the SM contribution corresponds to = 1 and all other couplings 
vanish, i.e. = Tsm- Of course, as in the SM case, each of the new contributions can 
be written as in eqs. (2.14) and (2.16)-(2.17), and their explicit expressions can be given 
in terms of the reduced form factors /, gi, g 2 and 1. Fortunately, it is not necessary to 
give the expressions of all of them, since they satisfy the following identities 


rpV 

-‘-LL 


rjnVQ 

RL 


rpVB 

-‘rl 


rpSV 

-‘lrrl 


rjiSV 

-‘llrr 


where {i = S,V,T) 


II- 

rpQ 

' LL1 

(3.4) 

A ''pSQ _ 
^-^RL - 

^ rpTQ ■ rpQ 

~ RL ~ RL^ 

(3.6) 

16T|f 

^rpTB ^rpST ■ rpB 

— g^RL ~ 2 RLRL — RL i 

(3.6) 

rpSV 

RLLR 

^rpVT '^rpVT ■ ^rpl 

~ Q LRRL ~ Q RLLR ~ ^ a i 

(3.7) 

rjiSV * 
-^RRLL 

- 2^/3, 

(3.8) 


Trl = rS + , (3.9) 


and analogous relations to (3.4)-(3.6) and (3.9) hold under L R. In order to read the 
results for these contributions eqs. (2.14) and (2.16) need to be used replacing by 
\Gui\‘^. 

In the limit of vanishing neutrino masses (i.e. keeping only left-handed neutrinos) the 
surviving structures are T^, and T^rrr- Consequently, gYi 7 ^ 1 and gf^R 7 ^ 0 would 
signal the most-likely first departures of the considered decays from the SM prediction. 

Fetscher et al. [17] have introduced a set of eight Michel-like parameters, Qij {i,j = 
L,R and YhijQij — 1)’ -Bhlj la and 1/3, which fully characterize the three-body 

decays of polarized muons in the case of the most general interactions dehned in (3.1). 
If we use this basis and the notation dehned in eqs. (3.4)-(3.8), then eq. (3.3) can be 
written as 


ng'-HP) 



+Se(/„ri + IfTf) 


(3.10) 


It is noteworthy that even for the polarized case, and keeping non-vanishing masses of hnal 
charged leptons, can be written in terms of the above eight Michel-like parameters. 
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We have checked that, in the limit of massless charged leptons in the final state, the last 
two terms in the previous expression vanish, in agreement with Eq. (4) in Ref. [18]. 

In the rest frame of the polarized decaying particle, and after integrating over neutrino 
four-momenta, the numerator in the r.h.s. of Eq. (2.13) has the general form (which inclu¬ 
des T-odd correlations) 




‘|G 




L Pi ■ S — Gi P2 ■ s — G2 P 3 • s 


+ HiS- (pi X P2) H2S- (p2 X P3) -F H^S- (pi X P3) 


(3.11) 


It is worth to note that the spin-independent coefficient F remains invariant under par¬ 
ity inversion, while the spin-dependent form factors L, Gi and G 2 (T-even correlation 
coefficients of spin and momenta) change sign under this operation. Correspondingly, F 
{L, Gi, G 2 ) has (have) the same (opposite sign) expression for the contributions and 
'^RR^ '^LLRR '^rrllj Therefore, we only need to provide five contributions, 

which we choose to be (SM case, see Appendix B), T^Yrr, and It 

must be noted that the latter two are a more convenient choice than and in 
order to get more compact expressions. However, both sets are related via 

rtf = (3.12) 

rff = i(i6r|i-rti). (3.13) 

A similar decomposition as in Eq. (3.10) can be written for the form factors in Eq. 
(3.11) ® 

F = {QllF?l + QrlFSl + BlrFEr + L^R)+ ^e{RFi + If,F^) , (3.14) 

L = {QllL^ll + QrlLIl + BlrLI^ + L ^ R) + , (3.15) 

Gi = {Q + QrlG^j^j^ + BlrG^i^p^ -|- L -h- i?) -|- -I- IpG\ p) ,(3.16) 

G 2 = + QrlG2,rl + BlrG2^i^r + A 0 R) -|- 'Sle[laG2,a + ^/sGi^js) ,(3.17) 


and 

Hi = H2 = H3 = 0 (3.18) 

even in the beyond the SM case. 

In Appendix D we provide the expressions for the reduced form factors /, /, gi and g 2 
(see their definition in terms of F, L, Gi and G 2 at the end of Section 2) for the 
T^Yrr, and structure. In the case of identical (non-identical) charged leptons, 
forty (twelve) different reduced form factors are required, in general, to specify each of 
the five independent Lorentz-chiral amplitudes defined above. As an illustration, we write 

®Or course, the symbol L O i? refers to the chiralities and not to the T-even correlation coefficient L. 
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the explicit expression of (analogous expressions hold for the others): 

4 

rriS \ ^ / nSAi iSAi -» -» SAi -» S,ii -* /o 1 r\\ 

^RL = 2^ Cii Pi ■ s - P 2 • S - P3 • Sj (3.19) 

2=1 

I \ ^ / rSAkki iSAkki -» -» SAkki -» -» SAkki -» 

+ ^ikCikki yjRL ~ ^RL Pi ' ^ — dlRL P2 ' S — (?2rl PS ' sj , 

i<k 

where cu, Cikki are the factors that relate the amplitudes Tjj, Tikki to the reduced ampli¬ 
tudes tii, tikki in Eq. (2.11), respectively, and Uik = +1 for {i,k) = (1,2), (3,4) and —1 
for {i,k) = (1, 3), (1,4), (2, 3), (2,4). Once all the independent amplitudes are obtained 
we can combine them using (3.4)-(3.9) and (3.12)-(3.13) and insert them into Eq. (3.10) 
to obtain the hnal expression for the squared amplitude. 

3.1 Possibility of T-odd correlation terms 

As it was pointed out before, the coefficients Hi (i=l,2,3) of the T-odd correlations may 
signal violation of time-reversal symmetry and are discussed in the following. T-odd 
correlations in the form of triple products of spin and momenta can be generated at 
higher orders due to exchange of virtual photons between charged leptons in five body 
decays [42] . These rescattering effects that mimic time reversal non-invariance would be 
suppressed by at least a factor of 0 {a/7r) with respect to the leading contributions already 
considered [42]. 

The contribution of the effective Lagrangian (3.1) to the coefficients Hi (i=l,2,3) van¬ 
ishes in the limit of massless neutrinos, due to the LL and RR structure for the product of 
weak currents in the amplitude (in the case of the interference between the tensor current 
and the 1/ ± A currents the result still holds but we could only check it by brute force). 
These coefficients would arise, in principle, from the imaginary parts of the product of 
couplings , i ^ j) in the interference terms of the squared amplitude; however 

they vanish for massless neutrinos, as already mentioned. Thus, interferences between dif¬ 
ferent currents giving rise to T violation will be extremely suppressed when considering 
(light) massive neutrinos for the underlying physics giving rise to the Lagrangian in (3.1) 
and we will therefore neglect it. 

We would like to emphasize that a measurement of T violation in these decays may 
not come from the Lagrangian (3.1). Particularly, if mixed chiral structures LR or RL 
for the product of weak currents are present, they would manifest as non-zero T-odd 
correlations in the energy and angular distributions. One may think, for instance, about 
the exchange of additional gauge bosons that couple to flavor-violating weak currents with 
both chiralities, as it would be the case for the contribution of a Z' neutral boson [43] 
(see also Sect. V.D of [44] and references therein) that couples non-diagonally to leptons 
as follows 

£ = ^ gz'fJiYi.Vij - + h-c.. (3.20) 

ij 
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In particular, it is possible that the interference of the corresponding amplitude with the 
SM V — A term could bring in T-odd terms. However, although allowing for the previous 
LFV vertices, the process L~ —)■ ulUi is nevertheless lepton flavor-conserving. 

This implies that Fierz rearrangements in the amplitude will bring the contribution of 
such interference into those corresponding to the structures already included into the 
lepton flavor conserving Lagrangian eq. (3.1). Consequently, the T violating form 
factors will again be zero, as we have checked by explicit computation. Precisely because 
of this general suppression of T-odd correlations in the SM and some of its extensions 
(by the tiny neutrino masses), searches for their effects in hve-body decays of leptons 
provide another place to look for (unexpected) indirect manifestations of CP violation 
in the lepton sector. On the other hand, their non-observation would be useful to place 
constraints (although probably mild) on the leptonic CP phases. 

At this point we shall recall that T violation has been studied in the LFV 

decays [45, 46] where these effects are not suppressed in general. Therefore, 
a non-vanishing measurement of T-odd correlations in leptonic tan decays with three 
charged leptons would most probably be an indirect signal of LFV neutrinoless processes. 

4 L — > 3£ 2z/ decays as backgrounds for LFV L ^ 3i searches 

The decays L~ are some of the main backgrounds in searches for LFV 

L~ decays as they may mimic the signal owing to undetected neutrinos. Our 

thorough analysis of the former processes allows to study how they would mask signals 
from the latter. In order to do this, we consider the study of L~ —)■ decays as 

reported in Ref. [47] (see Ref. [48] for an earlier effective-held theory analysis of these 
processes). There, a general low-energy ehective Lagrangian is considered to describe 
the LFV process and -in addition to the most popular dipole-type amplitude induced 
by radiative penguin diagram- it also incorporates ehective four-lepton operators with 
vector and scalar structures and the diherent allowed chiralities. Their equation (3.2) [47] 
gives the double diherential decay distribution in the two independent opposite-sign lepton 
invariant masses in terms of the coefficients of the relevant operators introduced above. 
From this, it is straightforward to obtain any observable of interest for our comparison. 

We have hrst reproduced Figure 5 of Ref. [47], for the same sign di-muon invari¬ 
ant mass spectrum in the r —?■ 3yU decay for which the following benchmark points are 
considered (see dehnitions in Ref. [47]): 

• Vector Model; Cvlr = Cvrl = 0.3 with all other couplings vanishing. 

• Scalar Model: Csll = Csrr = 1 with all other couplings vanishing. 

®The SM background to LFV r —>■ £7 processes given by the t —>■ i^vv decays is also computed in 

[46]. 

^°Their operators including two quark fields and their gluonic operators are not relevant for our dis¬ 
cussion, as opposed to LFV semileptonic tau decays. 
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• Dipole Model: Cdl = C'dj? = 0.1 with all other couplings vanishing. 

The scale of LFV is assumed to be A = 1 TeV. We have stick exactly to this setting 
in what follows. We have analyzed the same and opposite sign di-lepton invariant mass 
distribution and also the differential decay width versus E^- both for the SM L~ —)■ 
decays and for the L~ LFV processes. The best discriminating 

variable for this search of LFV turns out to be the opposite-sign lepton-pair invariant 
mass distribution. 

The comparison of the current upper limits on L~ —)■ decays to our SM 

predictions for the five-lepton lepton decays (see Table 1) suggests that the tan decay 
modes with two or three muons are better suited for these searches than those with two 
or three electrons. In the case of muon decays the huge background to signal ratio (i.e. 
the ratio among the SM process with neutrinos over the corresponding neutrinoless LFV 
decay) will make detection very challenging. 


Ghannel 

Gurrent upper limit (UL) [27, 49] 

S/B (UL) 

Expected UL [50] 

BR(r“ e“e+e“) 

1.4 • 10-® 

~ 3■10-^ 

~ 10-9 

BR(r“ —)■ e~ jji'^jji~) 

1.6 • 10“® 

~ 0.1 

~ 10-9 

BR(r“ —)■ 

1.1 • 10“® 

o 

1 

~ 10-9 

BR(r“ ^ 

1.2 • 10“® 

~ 0.1 

~ 10-9 

BR(/i“ e~e^e~) 

o 

o 

1 

to 

~ 3■10-® 

~ 10-16 


Table 2. Current and expected sensitivities on LFV L~ —>■ searches. The signal to 

background ratios (S/B) are estimated from current UL on BR’s of LFV decays (signal) and of 
five-body decays (background). 


In Fig. 4 we confront the SM background to the hypothetical signals of LFV according 
to the Vector, Scalar and Dipole Models in the case of three muons in the hnal state. The 
normalization of the new physics curves is chosen so that the branching ratios for signal 
over background give 0.1, which basically corresponds to a LFV signal about the current 
upper limit BR{t —)■ 3/i) < 1.2 ■ 10“® [49]. In this case, a cut for > 0.75 GeV^ 

will be most efficient. The similarity of the different new physics models above this cut 
suggests that it would be hard to disentangle the type of new physics with early data. 
If the new physics signal is set to the expected [50] upper limit ~ 10“®, the optimal 
cut moves to >1-2 GeV^. The r —)■ e/r/i case is very similar to the three muon 

channel, with a bit smaller cut: ^ 0.5 GeV^ for S/B 0.1 and ml+^. > 0.8 

GeV^ for the envisaged near-future upper limit. The cases with two or three electrons 
are much harder for detection due to the signal to background ratio < 10“^, which limits 
the region of the spectrum available for detection, as it can be seen in Figure 6 where a 
cut for > 1.25 GeV^ shall be needed (~ 1.5 GeV^ for r e/i/r). In addition to a 

very good statistics, an exquisite control of SM backgrounds will be needed to discover 
H 3e. 
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Figure 4. Comparison of opposite-sign lepton-pair invariant mass distribution in five-body and 
LFV decays of tau leptons. Three benchmark scenarios of new physics in r —)• 3/r are used, 
according to Ref. [47], and S/B ^ 0.1 (about the current upper limit) is assumed. 


5 Conclusions 

In this paper we have studied some aspects of all possible five-body leptonic channels of 
/i and r lepton decays. In all our calculations, we have kept the finite masses of the final 
state charged particles, but neglected neutrino masses. Firstly, we have re-calculated 
the branching ratios in the SM and have compared our results with two previous and 
conflicting calculations [20, 21]. Our results were obtained by using a different method 
for the integration of the five-body phase-space; our yields are in good agreement with 
the ones of Ref. [20] . 

In the second part, we have calculated the energy and angular distribution of the 
three-charged leptons, which can be written in terms of T-even and T-odd correlations 
of the spin s of the decaying particle and the momenta p* {i = 1,2,3) of final state 
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o 




Figure 5. SM background confronted to our benchmark scenarios of new physics in r —>■ 3e for 
S/B ^ 0.001 (about the current upper limit). 


charged leptons. We have derived the expressions for the T-even correlations of the 
form s • Pi, which are non-zero in the case of the SM and in the framework of the most 
general low-energy effective Lagrangian which conserves lepton-flavors. In particular, we 
have provided analytic expressions for the polarized decay probability in terms of the 
Michel-like parameters introduced by Fetscher et. al. [17] by keeping daughter-lepton 
mass dependence, as needed by Belle(-II). These results will be useful to test precisely 
the structure of the weak charged current in five-body leptonic decays of /i and r leptons. 

We have also shown that, for massless neutrinos, the coefficients of the T-odd cor¬ 
relations of the form s ■ (pj x pj) vanish in the SM and in the extension with the most 
general effective theory of the (current) x (current) form, even if we allow for lepton-flavor 
violating vertices with extra gauge bosons. Consequently, observation of non-vanishing 
T-odd correlations either will point to indirect signals of non-conventional mechanisms 
of CP violation in the lepton sector (if produced in five-lepton lepton decays) or to the 
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Figure 6. SM background confronted to our benchmark scenarios of new physics in /r —)• 3e for 
S/B ~ 0.001 (about the current upper limit). 


discovery of LFV (through the corresponding neutrinoless processes). These appealing 
characteristics, in our view, have not been noticed previously in the literature regarding 
purely leptonic interactions and open a new way to search for effects of CP violation. 

Finally, we have considered in detail the backgrounds that the SM L~ 
decays constitute to searches of LFV in the corresponding neutrinoless decays. Tan decays 
with two or three muons are favoured with a moderate signal to background ratio. Gen¬ 
erally, the differential decay distribution versus the invariant mass of the opposite-sign 
lepton pair is the best observable for maximizing the signal region. This looks similar 
to the discriminating power of the di-lepton energy spectrum in searches of neutrinoless 
nuclear double beta decays which would show up as a single line excess over the con¬ 
tinuum distribution due to the allowed two neutrinos double beta decays. Forthcoming 
experiments will be able to perform these searches with reasonable cuts on this variable 
in the sub-GeV region. 
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Appendix A 

We include here the expressions for the reduced amplitudes in the Standard Model case 
for unpolarized tau decays, as defined in Eq. (2.11). For simplicity in what follows we 
will use the following notation, Qi = Q ■ Pi and pij = pi ■ pj [i ^ j). All results shown in 
the Appendices were obtained using FeynGalc [51]. 

til = - 2 {ml + Pi2 + Pl 3 + PlA - Ql) {P 23 P 2 aQ 2 + PsaQI + P 23 P 3 aQ 3 “ P2AQ2Q3 ( 5 . 1 ) 
-P3AQ2Q3 + P2aQI + + P23) {p2A + P 3 A - Q4) - TTI^Qa - ( - 2Q2Q3 + P23 

X (Q2 + Q 3 )')Qa + {P 2 a{ 2 Q 2 + Q3) + P 3 a{Q 2 + — (^23 + Q2 + ^3)^4)) ; 

t22 = 2 (mi(m^ P23) -h m^( - 2 pi 2 - 2 pi 3 -F p^ - Qi) - ^13(^13(^23 - P2a + Q2) 

+P 23 {PlA + P23 + P 3 A -Q1-Q3)) - Pl 2 ( 2 pi 3 ^23 (Pl 4 + P 23 + P2A - Ql - Q2) 
+Pl 3 ( 2 pi 4 + 4 p 23 + P 2 A + P34 ~ 2 Qi — Q2 — Q3)) + ml ( 3 m^ + P12 ( - 2pi3 -F P23) 
+P23 {Pl 3 + Pl 4 + 2 p 23 + P 2 A + P34 — Ql — Q2 — Q3) + {Pl 2 + Pl 3 + Pl 4 + 5^23 
+P24 + P 34 - Ql - Q2 - Q3)) - P? 2 ( 2 P 13 + P23 - P 3 A + Q3) 2 pl 2 - 2 p\^ 

+P 23 {P 1 A - Ql) - Pl 2 ( 6 pi 3 + Pl 4 + 3 p 23 + 2^24 + P 3 A - Ql - 2(^2 “ Q3) 

+Pl3( ~ Pl4 ~ 3p23 ~ P 24 ~ 2p34 + Ql + Q 2 + 2 Q 2 ,)) )Q4 , (5-2) 

tl221 = —2 (Pi3 ( — P34Q2 + P24Q3 + 2Q2Q4) + P12 {P 3 aQ 2 — P24Q3 + 2Q3Q4) -|- ml ( (pi2 
—P13) {P 3 aQ 2 — P24Q3) + ~ PlA {Q2 + Q3) + Ql {P 2 A + P34 — 2Q4) ) -|- P23 

X (Ql (P 24 + P 3 A — 2Q4) — Pi4 [Q2 + Q3) ) + 2 (P13Q2 + P12Q3) Q4) + (Ql (P 24 
+P 34 ~ 4 Q 4 ) — Pi4(Q2 + Q 3 + 2 Q 4 )) -|- M^(2m^Pi4 -|- P23(P13P24 + P12P3A — {Pl 2 
+P13) Qa) + { 2 piAP 23 - P12 {p2A + Qa) - Pl 3 (P 34 + Q4) ) ) + Pl 3 ( (Pl 4 + 2^24 
—Ql — 2Q2) ( — P34Q2 + P24Q3) + (( 2 pi 4 + P24 ~ 2 Qi — Q2) Q2 + ( ~ P2A + Q2) 
xQ 3 )Q 4 + P 23 (P 34 (Qi ~ Q2) ~ P 14 (Q 2 + Q3) + P 24 (Qi + Q3 ~ Q4) + Q 4 ( ~ 2 Qi 
-I-Q2 -l- Q4))) + P12 ((pi 4 + 2p34 — Ql — 2Q3) (P34Q2 — P24Q3) + ( ( 2 pi 4 — 2 Qi 
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+^2 ~ Qs) Qs + P34 { — Q2 + Q3) + 2pi3 (Q2 + Qs) ) Q4 + P23 (p24 {Qi — Q3) 

—pu {Q2 + Q3) + P34 {Qi + Q2 — Qa) + Q4 ( ~ 2Qi + Qs + Q4) ) ) — P23 {Pl 4 {Q2 
+Q3) + Pl4( ~ Q1Q2 + P34( ~ Ql + Q2) — Q1Q3 ~ 2Q2Q3 + P24( — Ql Q3) 
+^23(^2 + Q3) + 2Q1Q4) + Ql(p 34(^1 + Q2) — P 23 {P 24 + P34 ~ 2Q4) + (p34 

— 2 Qi — Q2 — Q 3 )Q 4 +P 24 ( ~ 2^34 + Ql + Q3 + Q4))) + '1^'^ {P 13 P 24 :Qi 
+P14P24Q1 + 2P23P24Q1 - PI4Q1 + P13P34Q1 + PIAP34Q1 + 2P23P34Q1 - PuQl 
-P 24 Q\ - P 3 iQ\ - P13PUQ2 - PI4Q2 - 2P14P23Q2 + P14P24Q2 - P13P34Q2 
+P14Q1Q2 + P24Q1Q2 — P14Q2 ~ P13P14Q3 ~ P14Q3 ~ 2P14P23Q3 + P13P24Q3 

+P14P34Q3 + P14Q1Q3 + P34Q1Q3 - P14Q3 + ( - 2pi4(p23 + Ql) + Ql{- GP 23 

—P 24 — P34 + 2Qi + Q2 + Q3) + Pl3 (p34 ~ 2Qi + 2Q2 + Qs) ) Q4 + P13Q4 
+P12 (p34 (Ql + Q2) ~ Pl4 {Q2 + Q3) + P24 {Ql — Qs + Q4) + Q4 ( ~ 2Qi + Q2 
+2Q3 + Q4)))), ( 5 - 3 ) 

^1331 = 3m® Q 4 — m'^(pi 4 ( 2 Qi + Q 2 + 2 Q 3 ) + 2^34 (Qi + Q 2 + 4 Q 3 ) + p24(Qi + 2 (Q 2 
+Q3)) ~ (pi2 + 3pi3 + 3p23 ~ P34 + 5Q3)Q4 — 2Q4) + 2Q3(pi2P34 + 2p23P34Ql 
+P24P34Q1 ~ 2P34Q1Q2 ~ P12P34Q3 + P24Q1Q3 ~ ^13(^34 (2p23 + P24 ~ 2Q2) 
+P24Q3) + Pl4 ( — 2P24Q3 ~ P23 {P 34 + Q3 — Qa) + Q2 (P34 + Qs — Q4) ) + Pl3 

X (2^23 + P24 — 2Q2)Q4 + ( ~ Ql (2p23 + P24 — 2Q2) + Pl2( — 2^34 + Q3))Q4 

+P12Q4) ~ {P23P34Q1 — P24P34Q1 — P34Q1 + P24Q1 + P34Q1 + 2P13P24Q2 

+P13P34Q2 - PI4Q2 - P24Q1Q2 + P34Q2 + P13P24Q3 + 2P12P34Q3 + GP13P34Q3 
+6P23P34Q3 + P24P34Q3 + P24Q1Q3 ~ 4P34Q1Q3 — P24Q2Q3 — 4P34Q2Q3 + 4P24Q3 
+(^23(^34 — Q2 — 6Q3) + ^12(^34 ~ Ql ~ Q2 ~ 3Q3) — ^13(2^23 ~ P34 + Ql 
+6Q3) + 2 (p 24 Qi + P34(Qi + Q2) + Q3(3(Qi + Q2) + Q3)))Q4 ~ 2(pi2 + P13 
+P23 ~ Q 3 )Q 4 + Pl 4 ( ~ P34Q2 + Q2 + P34Q3 + Q2Q3 + 4Q3 + P 23 ( 2 Qi + Q3) 

~Ql (Q2 + Q3) + 2Q2Q4) ) + (2^34 ( — P14P23 + P12P34 + P23Q1 + Pl3 ( ~ 2^23 

-P24 + Q2)) + (P14P23 -P23Q1 +Pl3(2p23 + P24 “ Q2) “ Pl2 (P34 + Q3))Q4 

+m'^(pi4 + P24 - 2p34 + 3Q4) + m^ (P14P23 + 2pi4P24 - P 23 P 2 A “ 4^23^34 + 2p34Ql 

+2P34Q2 ~ Pl 2 (Pl 4 + P24 + 2^34 — Q4) + (pi4 + 3 p 23 + P 24 + P 34 ~ Ql ~ Q2 
+Q3) Q4 + Pl3 ( - Pl4 + P24 - 4p34 + 3Q4) ) ) , (5.4) 

^1441 = {m^ (2pi4 - P24 + P34) + {P 23 (2pi4 - P24 + P34) - Pl3 (P24 + P34) + Pl2 

X ( - P24 + P34 - 2Q4)) + 2 pi 2 P 23 {p 3 A “ Q4)) + m^(3p34Ql “ P14Q2 “ 8P34Q2 

— 3P14Q3 — (2pi4 +P34 + 6Q1 — 9Q2 + 9Q3)Q4 +P24 (Qi + 8Q3 + Q4)) + 

X (2P12P34Q1 + 2P13P34Q1 + 2P14P34Q1 - P^4 Qi - 2p34Q^ - 5P12P34Q2 - 9pi3P34 
XQ2 — 3P14P34Q2 — P34Q2 + 3P34Q1Q2 — P14Q2 + P34Q2 ~ 2P12P14Q3 ~ 2P13P14Q3 
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—2pl^Q^ + PiiPsiQs + ‘^PiiQiQs + PsiQiQs + P34Q2Q3 ~ PuQl + P24{ ~ Qi 
+Q 3 ) ~ ( “^ 13(^34 ~ 2(5i + IIQ 2 ~ 5 Q 3 ) +^ 12(^34 + 2Qi — 7 Q 2 + 3 Q 3 ) 
+2 (pi 4 Qi — Qi — 2P14Q2 — P34Q2 + Q1Q2 + Q2 + (2pi4 + P34 ~ 2 Qi)Q 3 — Q3)) 
XQ4 + 2pi2Q‘l + P24( ~ P34Q2 + Q1Q2 + 5P12Q3 + 9^13^3 + P34Q3 ~ 3Q1Q3 
— Q2Q3 ~ Q3+ Pl 4 {Q2 + 3Q3) + (pi2 + Pl 3 ~ 2 (Ql — Q2 + Q3) ) Q4) + P23 (5^34 
xQl- P14Q2 - 4^34^2 - SPmQs - (2pi4 + P34 + 8Q1 - 5Q2 + 5 Q 3 )Q 4 + P 24 {Qi 
+4Q3 + Q^j ) ) + 2 [PI2Q3Q4 + {pi 3 + Pl 4 + P23 + P24 — Ql — Q2) {p 23 {P34Q1 
—P14Q3 ~ Q1Q4) + Pl3 ( ~ P34Q2 + P24Q3 + Q2Q4) ) + P12 ( ~ {Pl 3 ~ P34 + Q3) 

X {P34Q2 - P24Q3) + ( (Pl3 - P34) Q2 + (Pl3 + Pl4 -Ql+ Q2) Qs) Q4 + P23 {P34Q1 

—PuQs ~ {P 34 + Qi)Q 4 + QT))) ) ( 5 - 5 ) 

^2442 = -4 (2m® + m^( - 4pi2 + 3pi3 + pu + 3 p 23 + P24 - Qi - Q2) + 3pi2 

+ {pi3 + P23) [plS + Pl4 + P23 + P24 - Ql- Q2) - 2,pi2 (2pi3 + Pi4 + 2p23 + ^24 
+P34 — Ql — Q2 — Q3) ) ~ Pl2 ( (pi 3 + P23) {pi3 + Pl 4 + P23 + P24 — Ql — Q2) 


+P 12 {Pl3 + P 23 - P34 + Q 3 ) ) ) Q 4 , 


(5.6) 

1 , 

^44 — ^29 1 

^z|pi^p2,7ni^7n ’ 


t2332 =tl44l\p^^p^, 

tT44T — tl221 

(5.7) 


Appendix B 

Here, we collect the expressions for the reduced spin-(in)dependent form factors in the 
Standard Model case for polarized tans (we recall that the SM contribution, Tsm, equals 
the contribution of the most general low-energy effective held theory). The origin of 
the identities below can be traced back to the exchange of indistinguishable fermions as 
well as to the symmetric character of the electromagnetic tensor p^pj -|- pjp'^ — g^^{pi ■ 
Pj + m?), which comes from summing over the polarizations of the pair. 

fl]^ = 2{2{ml + pi2 + pi3 -Qi){{m'^ + P23 - Q2) {M‘^P23 + + Q2)) 

+ ((m^ -|- P23) ~ 2 (m^ -|- P23) Q2 + 2Q2) Q3 + + 2 (ml + pi2 + P13 — Qi) 

X (M^ (m^ P23) + (m"^ -h m^p23 - 2Q2Q3) - P23 {Ql + Ql) - 2 m^ {Ql 

+Q2Q3 + Q3)) + 2(m^ -|- P12 + P13 — Ql) {M‘^{m? + P23) (m^ -|- P23 — Q3) 

-|-m'^ {Q2 + Q3) + Q2 {P23 — 2P23Q3 + 2Q3) -|- m^ {Ql — 2Q2Q3 — Q3 + P23 (2Q2 
+Q3))) + (ilT^ -|- 3 ml + 2 {m? + 2pi2 + 2pi3 + P23 — 2 Q\ — Q2 — Q3)) {M‘^{m?‘ 
+P23) {pi2 + P13 — Ql) — m^Qi -|- P12P23Q2 — P23Q1Q2 + P13QI + P13P23Q3 
—P23Q1Q3 — P12Q2Q3 — P13Q2Q3 + 2Q1Q2Q3 + P12Q3 + { — Ql {p 23 + Q2 

+Q3) + P12 {‘2Q2 + Qs) + Pi 3 {Q2 + 2Q3)))) , ( 5 . 8 ) 
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~ 2 iVf ( 2 m^ + M‘^ + + 4 pi 2 + + 2 p 23 ~ 4 Qi — 2Q2 ~ ‘^Qs) 

X + P23) ~ 2(52(53) ) ( 5 - 9 ) 

5 'ism = -2M(m'^(8pi2 + 6^13 - QQi) + M^(pi 2 P 23 + m^(2pi2 +P13 - Qi) - P23Q1 
+Pi 3 iQ 2 -Qs)) + + 3 pi 2 P 23 + 2^23 + ”^^(6^12 + 3 pi 3 + 6^23 - 3 ( 5 i 

—4(52 — QQs ) + 3pi3((52 — Qs) + 2(53(^52 + (53) ~P23(3(5 i + 2(52 + 4(53)) 
+ 2 m^( 4 p ^2 + 2 Pi 3 + Pl 2 ( 6 pi 3 + 6 p 23 - QQl - 4(52 - 5 ( 53 ) + Pl 3 ( 4 p 23 - 4 ( 5 l 
— 2Q2 — 5(53) + ( 5 l( ~ 5^23 + 2 ( 5 l + 3(52 + 4(53)) + 2 ( 2 p^ 2 P 23 + 2 pi 3((52 ~ ^3) 
+Pl 3 ((p 23 ~ 2 ( 5 l — Q2) {p 23 + Q2) + ( ~ 3^23 + 2 ( 5 l + Q 2 )Q 3 + 2 ( 5 ^) + Ql 
2^23 + 2^23 (Ql + Q2) + 3^23(53 ~ <53 (<52 + Q3)) +^12(2^23 + 2^13(^23 
+(52 ~ (53) + Q 3 {Q 2 + (53) ~ P 23 ( 4 ( 5 i + 2(52 + 3(53)))) ) ( 5 . 10 ) 

fsM ~ — 2 ( 2 m^Qi + m^( 3 m^(( 5 i + <52 + ^3) + 3 ^ 23 (Qi + Q2 + Q3)) ~ 2 m‘^( 2 Qi 

X ( — P23 + Ql + <52 + Q3) + P12 (2Q1 + 4Q2 + 3Q3) + P13 (2Q1 + 3Q2 + 4Q3)) 
— 2 m^ (pi2 (3Q1 + 4Q2 + 2Q3) + p\^ (3Q1 + 2Q2 + 4Q3) + P23Q1 ( ~ P23 + 2 (Qi 
+Q 2 + Q3)) + 2^13(2^23(Ql + Q2) + 3P23Q3 ~ (Ql + Q2 + Q3) (Ql + Q2 
+2Q3)) + 2 pi 2 ( 5 pi 3 (Ql + Q2 + Q3) ~ (Ql + Q2 + Q3) (Ql + 2Q2 + Q3) + P23 
X (2Q1 + 3Q2 + 2Q3))) + rn\ (P13Q2 + P12Q3 + P13P23 (3Q1 + 4Q2 + 3Q3) 
+m^( 9 Qi + 8(Q2 + Q3)) + 2 p 23 (Qi + Q2 + Q3) ( 3 p 23 ~ 2 (Qi + Q2 + Q3)) 
+m^( 3 pi 3 Qi + I5P23Q1 ~ 4 Q^ + 3P13Q2 + I4P23Q2 ~ 8Q1Q2 — 4Q2 + 2 (pi 3 
+ 7^23 ~ 4 (Qi + Q2))Q3 ~ 4Q3 + Pi2(3Qi + 2Q2 + 3Q3)) + Pl 2 ( 3 p 23 (Ql + Q2) 
+4P23Q3 - Pi 3 ( 6 Qi + 5 (Q 2 + Q3)))) + - 2 m^ - 2^23) + ’^^( 4 (pi 2 

+P13) + 3 Qi) + ^12(4^13 + 2^23 + 3Q3) + ^13(^13(2^23 + 3Q2) + ^23(2^23 

- 3 (Qi + Q3))) + ”^?( ~ + P12( 4 pi 3 - 2^23) + m^( - 2pi2 - 2pi^ - 10 p 23 

+3 (Ql + Q2 + Q3)) + P23 ( ~ 2pi3 — 4 p 23 + 3 (Ql + Q2 + Q3))) + P12 ( 4 pi 3 
+P23 ( 2^23 ~ 3 (Qi+Q 2 ))+ Pi 3 ( 8^23 ~ 3 ( 2 Qi + Q2 + Q3) ) ) + ( 4 pi 2 + 4 pi 3 

+3P23Q1 + 3^12(4^13 + 2 p 23 — Ql — 2Q2 — Q3) + ^13(6^23 ~ 3 (Qi + Q2 + 2 
XQ3)))) + 2 (Pi 2 Q 3 — Pi 2 (P 13 ( 4 Qi + 3Q2 + 2Q3) + 2 (p 23 (Qi + Q2) + Q3 (Ql 
+Q2 + Q3) ) ) + Pl 3 (P13Q2 ~ 2pi3 (p 23 (Ql + Q3) + Q2 (Ql + Q2 + Q3) ) + P23 (2 
X (Qi + Qs) (Qi + Q2 + Q3) ~ P23 (2Q1 + Q2 + 2Q3))) + P12 ( ~ Pi 3 (4Q1 + 2Q2 
+3Q3) + P 23 ( 2 (Qi + Q2) (Ql + Q2 + Q3) ~ P 23 ( 2 (Qi + Q2) + Q3)) + 2 pi 3 ((Qi 
+Q2 + Q3) (2Q1 + Q2 + Q3) ~ P 23 ( 4 Qi + 3 (Q 2 + Q3)))))) ) ( 5 - 11 ) 

= - 2 M ( - 2 m® - 3 mf (m^ + P23) + ( - m'^ + 2 pi 2 Pi 3 - mj (m^ + P23) 

+m^(pi2 + P13 ~ P23) + (pi2 + ^13)^23) + 2 m^( 2 pi 2 + 2 pi 3 — 2 p 23 + 2 Qi + Q2 
+Q3) + 2 m^( 3 p ^2 + 3 pi 3 + ^13(4^23 ~ 2 (Qi + Q2) — 3Q3) + pi 2 (l 0 pi 3 + 4 p 23 
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— 2 Qi — 3Q2 — 2Q3) + P23 ( — P23 + 2 Qi + Q2 + Q3 )) + ~ 9 m^ + pi2 ( 6 pi 3 

~ 3 p 23 ) + — 3 pi 2 — 3 pi 3 — 15 p 23 + 4 (Qi + Q2 + Qs)) + P 23 { ~ 3 pi 3 — 6 p 23 

+4 {Qi + Q2 + Q3))) +2 (pi 3 (pi 3 ( 2^23 + (52)+ P23 ( 2 p 23 ~ 2 ( 5 l — Q2 — 2Q3) ) 
+P12 ( 4 pi 3 + 2 p 23 + Q3) + P12 ( 4 pi 3 + P23 ( 2^23 ~ 2 (Qi + Q2) — Q3) + Pl 3 {^P 23 
— 4 Qi — 3(Q2 + Q3))))) , ( 3 - 12 ) 

gifj^ = 2 M {mj ( 3 m^ + 8^23) + (2pi2 + pn) + Pn + mj (m^ + P23) - P12 {pi 3 

+P23)) - 2 m‘^[Api 2 + 3 pi 3 + Qi) +m‘^{- ^{ 2 pl^ + 5pi2Pi3 + PI3 + ^Pi 2 P 23 
+2P13P23) + 2 ( 3 pi 2 + 2 pi 3 - P23) Qi + 4 [ 2 pi 2 + P13) Q2 + 6 (pi2 + P13) Q3) 

+m\ ( 8 m'^ + Pi3 + 4 pi 3 P 23 + P12 ( - ^Pi 3 + 3^23) + rn^ {‘2pi2 + ^Pis + 14^23 
—4 (Qi + Q2 + Qs) ) + 2^23 ( 3^23 — 2 (Qi + Q2 + Q3) ) ) + 2 ( — P12 ( 3 pi 3 + 2^23 
+Q3) + Pl 3 {Pl 3 + P23 ( ~ P23 + Ql + Q3) ~ Pl 3 {Qi + 2Q2 + Q3) ) + P12 ( — 2^33 
+P 23 ( 2 (Qi + Q2) + Qs ~ 2^23) + Pi 3 { ~ 6 p 23 + 3 Qi + 2(Q2 + Q3))))) ) ( 5 . 13 ) 

IFm = 2(M'^(3m^(pi2+P13) + 3 (pi 2 +^13)^23) + 12 m®( 5 i + ( - 3m^(Q2 + Qs) 

~ 3 p 23 {Q2 + Q3)) + 4 m^ (pi3 (4Q1 — Q2) + P12 (4Q1 — Qs) — 2(^1 ( — 4 p 23 + 3 Qi 
+2(Q2 + Qs))) + 4 ? 7 i^(pi 3 (Qi — 2Q2) + ^12(^1 ~ 2Q3) + ^13(7^23^1 ~ 3 Qi 
—2P23Q2 + Q1Q2 + 2Q2 ~ Q1Q3) + P12 ( ~ Qi ( ~ 2pi3 — 7 p 23 + 3 Qi + Q2) 

+ (Qi ~ 2^23) Q3 + 2Q3) + Ql ( 7^23 + {Qi + Q2 + Q3) {2Q1 + Q2 + Qs) ~ P23 
X (lOQi + 7 (Q2 + Qs)))) + 4 ( — P13Q2 + P23Q1 ( 2 p 23 ~ 2 ,Qi — Q2 — Q3) {P23 
— Ql — Q2 — Q3) ~ P12Q3 + Pl3 {p 23 {Ql — 2Q2) + Q2 (3 {Ql + Q2) + Q3) ) + P12 
X (P13Q2 + P23 {Ql — 2(^3) + Q3 (3Q1 + Q2 + 3(^3) ) + Pl 3 {p23 (3Q1 — Q2) — Q2 
X (2Q1 + Q2 — Qs) {Qi + Q2 + Qs) + P23 (2Q2 ~ ^Ql ~ ^QlQs) ) + P12 (p 23 (3Q1 
“Qs) + P13Q3 ~ Q 3 ( 2 Qi — Q2 + Q3) {Ql + Q2 + Q3) + P 23 ( ~ ^Qi ~ 2Q1Q2 
+2Q3) + Pl 3 ( 2 p 23 Ql + (Q2 ~ Q3) + Qi(Q2 + Q3)))) + M‘^{Am^{pi2 + Pl 3 

+ 2 Ql) +Pi 3 ( 2 p 23 ~ 5Q2) +^12(2^23 ~ 5Q3) + 3 ^ 23^1 ( 2^23 ~ 2 Qi — Q2 — Qs) 

+m\ {Am^ + P12P23 + P13P23 + 4 p 23 Qi - ^Pi 3 Q 2 - P23Q2 + {pi2 + P13 + 4 p 23 
+ 4 Qi — Q2 — Q3) — ( 4 pi 2 + P23) (53)+ P12 ( 4^23 + P23 {Ql — 4Q2 — 8Q3) + 3Q3 
X {2Q1 — Q2 + Qs) + Pl 3 ( 8 p 23 ~ 3 {Q2 + Q3) ) ) + {pi2 ( 4 pi 3 + 8 p 23 + Ql — 4 Q 2 
— 8Q3) + Pi3 { 8 p 23 + Ql — 8Q2 — 4Q3) — Ql{ — 14^23 + 3 {2Q1 + Q2 + Qs) ) ) 

+^13(4^23 + 3 Q 2 ( 2 Qi + Q2 — Qs) + P 23 {Qi ~ 4(2^2 + Qs)))) + ^l( ~ 3P33Q2 

—8PI2Q3 + ^^( 3 pi 2 Qi + 3 pi 3 Qi + 14 ^ 23^1 ~ SQi — 4P12Q2 ~ 8P13Q2 ~ 8P23Q2 
—Q1Q2 — {8pi2 + 4 pi 3 + 8 p 23 + Ql — 4 Q 2 )Q 3 ) + Tn‘^{8Qi — 4(^2 + Qs)) + P23 
X ( — 6 Q^ — Ql (Q2 + Qs) + 2 (Q 2 + 4Q2Q3 + Qs) + P 23 ( 6 Qi — 4 :{Q 2 + Qs))) 
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+P 13 (Q 2 ( 6 Qi + 5Q2 + 3Q3) + P 23 ( 3 Qi — 4 ( 2 Q 2 + Q3))) + Pl 2 ( 3 pi 3 (Q 2 + Q3) 
+Q3(6 Qi + 3Q2 + 5(53) + P 23(3 Qi — 4(Q2 + 2Q3))))) ) (5-14) 

(l 2 m® + - P12 - P13 + 10^23 - 2(5i) - ^23(^13 - 2 p 23 

+ 2 Qi) + 2 P 13 Q 2 — Pi2{p23 ~ 2 Q 3 )) + 4m'^(4pi2 + 4pi3 + 8p23 — 3(2Qi + Q 2 

+Q3)) + 2771^(2^^2 + 2pi3 + 14^23 ~ 2OP23Q1 + 4 :Ql ~ IOP23Q2 + 4Q1Q2 ~ Q2 
+^ 13 ( 14^23 — 6 Q 1 + 2 Q 2 — Qs) — IOP 23 Q 3 + 4:QiQ3 — QI+ ^ 12 ( 4^13 + 14p23 
— 6 Q 1 — Q 2 + 2 Q 3 )) + ml( 8 m^ + 3 ^ 13(^23 + 2 Q 2 ) + 2 ^ 23 ( 3^23 ~ 3Qi — Q 2 
— Q3) + 3^12(^23 + 2Q3) + 7?7^(3pi2 + 3pi3 + 14^23 ~ 2(3Qi + Q2 + Q3))) 
+2(2p^3(p23 + 2Q2) + 2^32(^23 + 2Q3) + 2^23(2^23 + 2(^1 (Qi + Q2) + (2Q1 
+Q2) Q3 ~ 2^23 (2Q1 + Q2 + Qs) ) + Pl3 (6^23 + P23 ( ~ SQl + Q2 — ‘2Q3) 

— Q 2 { 4 :Qi + 3Q2 + Qs)) + ^12(6^23 + P 23 ( ~ QQl ~ 2Q2 + Qs) + 4^13(^23 + Q2 

+Q3) ~ Q3(4 Qi + Q2 + 3Q3)))) , ( 5.15) 

gi]?M = ‘2M{m\{8m^ + 3^23) + 4 m'^(pi 3 + Qi) + + 8^13(^23 - Q2) + 2Qi(pi2 

+ 4 p 23 - 2 ( 5 l - 8Q2 - 2(53)) + 7713(4777'^ - 3 pi 2 Pl 3 + 3^13 + 4 pi 2 P 23 + ^Pl 3 P 23 
+4^23 ~ P 23 Q 1 ~ 5 P 13 Q 2 ~ 2 P 23 Q 2 + 7?7^(4pi2 + 8pi3 + 8^23 ~ Ql ~ 

- (2pi2 + Pi 3 + 4^23) Q3) + ( - P12 (pi 3 - 2^23) + rnl [m^ + P23) + P23Q1 

+777^(4pi2 + 2pi3 + Qi) + Pl3(pi3 — Q 2 + Q 3 )) + 2 ( 2^33 + P 13 (4p23 ~ 2 ,Qi 

—6Q2 — Qs) + P23Q1 (2p23 ~ 2 (Qi + Q2) — Q3) ~ P12 (2pi3 + Q3) + P12 (2P23Q1 
+ {Qi — Q2 — Q3) Q3 + 2pi3 (Qi — Q2 + Q3) ) + Pl3 (2^23 + P ‘23 {Qi ~ 4Q2) 
+3Q1Q2 + (2Q2 ~ Q3) {Q2 + Qs)))) ; (5.16) 

5 ' 25 M = 9 is^m\p 2 ^p^ ; fi'asM = dl'^Allp^^p^ 5 92 ^sM = S'iFm Ip 2 «.p 3 ’ (^• 1 '^) 
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9ism 

pi<r^P2,Tni^m ’ 

92 SM 

y^SM \p^^p2_^rni^7n ’ 
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i44 
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^ Q 22 1 

y ^SM \p-^^p2^rni^m ’ 

(5.19) 
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— ‘sm| 


44 
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9ism 


92 SM 

y^SM \p^^p2_^rni^7n ’ 



/sM^ ~ 7 ? 7 ®( 13 (Qi + Q2) + 45Q3) + 7 ? 7 '^(l 0 pi 2 Ql + 8pisQl + 2IP23Q1 ~ QQl + IOP12Q2 

+2IP13Q2 + 8P23Q2 — 32Q1Q2 — QQl + ( 25 pi 2 + 64 pi 3 + 64^23 ~ 68(Qi + Q2)) 
xQs — 28 QI) — 8Q3( — 2(pi3 — Qi) {j)2s — Q2) {pis + P23 — Qi — Q2 — Qs) 
+Pl2( — (Pl3 ~ Ql) {P 23 ~ Q2) + ( ~ Pl3 ~ P23 + Ql+ Q 2 )Q 3 + Qs)) + 

X (12777'^ + 8 pi 3 P 23 - 3 p 23 Ql - 8P13Q2 - SpuQs + 777 ^( 4 pi 2 + 3 ( 4 pi 3 + 4^23 “ Ql 

— Q 2 + Q 3 ) ) ) + 2777^ (Ql (4^23 + 4 Q 2 (Qi + Q 2 ) ~ 3]723 (Qi + 3 Q 2 ) ) ~ PuQs 
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+ 12(^23 ~ 3^23^1 + Ql ~ 2^23^2 + 3QiQ2 + ^2)^8 + 12( — P23 + Ql + <52)^3 
+2Q3 + 4p^3(Q2 + SQs) +^12(5^23^1 + 5P13Q2 — 8Q1Q2 + (8pi3 + 8p23 — 9((5i 
+Q2))Q3 + 6Q3) + ^13(4^23(^1 + Q2) + 34^23^3 ~ 3(3QiQ2 + Q2 + 8Q1Q3 
+ I 2 Q 2 Q 3 + 4 Q 3 ))) + iVf^(24m^ + m^(24pi2 + 46pi3 + 46p23 ~ 20 (Qi + Q 2 )) 
+2m^(4p32 + llPl3 + 30pi3P23 + llpL ~ 12pi3(5l — 13^23^1 + 2Qi — I 3 P 13 Q 2 

— 12^23^2 + 2Q1Q2 + 2Q2 +^12(11^13 + llp23 ~ 4 (Qi + Q2) — 3Q3) — (3pi3 
+3p23 + 7 {Qi + Q2))Q3 ~ 4Q3) — 2(p^3( — 8p23 + 3Q2) + P12Q3 +Pl3( ~ 8^23 

+9^23 (Ql + Q2) — 2Q2 (Qi + Q2 ~ Qs) + 4P23Q3) + Ql(3p23 ~ 2^23 (Ql + Q2 
“Qs) ~ 6Q2Q3) + P12 (P23 (Qi + Qs) + Pl3 ( ~ 4p23 + Q2 + Q3) ~ 2Q3 (Qi + Q2 
+2Q3)))), (5.20) 

= M{- 13m® + M^(m^ + P23) + m^( - 10pi2 - 8^13 - 21^23 + 6Q1 + I6Q2 
+ I8Q3) +4Q3(p23(Pl2 + 2(pi3+P23-Ql)) “ (Pl2 + 3 (pi3 + ^23 “ Ql) ) Q2 
+Q2 + (pi2 ~ 2^23 + 2Q2)Q3) ~ 2m^ (5^12^23 + 4^13^23 + 4^23 ~ 3p23Ql 
—4P12Q2 — 2P13Q2 — 7P23Q2 + Q1Q2 + 3Q2 + ( ~ 3pi2 — 5pi3 — 13p23 + 5Qi 

+ 12Q2)Q3 + 6Q3) + M^(4m^ + 2^23(^12 + 3pi3 +P23 ~ 2Qi) — 2(pi3 +P23)Q2 

—2(pi2 + Q 2)Q3 + 2m^(pi2 + 3pi3 + 3p23 ~ 2Qi — Q2 + Q3))) > (5-21) 

fl'25M = M{- 45m® - MVi2 - 2M ^ {pl^ + 4^13^23 - 2P23Q1 + P12 (pi3 + P23 - Qi - Q2) 
~2pi3Q2) + m"^( — \2M‘^ — 25pi2 — 64pi3 — 64^23 + 50 (Qi + Q2) + 28Q3) 
+m^{M^ - 2 M^{pi 2 + 5pi3 + 5p23 - 2(Qi + Q2)) + 2(pi2 - 8^12^13 - 12pi3 
—8P12P23 ~ 34 pi 3 P 23 ~ 12^23 + 6P12Q1 + lOpisQi + 23P23Q1 ~ 7 Q'l + QP12Q2 
+23P13Q2 + I9P23Q2 ~ I2Q1Q2 — 7Q2 — 6(pi 2 — 2pi3 — 2p23 + Ql + Q2)Q3 

— 2Q3)) — 4(2(2^13^23 — P23Q1 — P13Q2) (pi3 + P23 — Ql — Q2 — Q3) — ^12(^23 

^ (Ql ~ 2Q3) + Pi3( — 2p23 + Q2 — 2Q3) + Q3(Qi + Q2 + 2Q3)))) , ( 5.22) 

fsM^ = -2(M'^(3mV2 + 3P12P23) + m®(lOQi - I5Q2 + 8Q3) + m^(7pi2Qi + 12pi3Qi 
+23P23Q1 ~ 15Qi — IOP12Q2 ~ 3IP13Q2 — 23P23Q2 + IOQ1Q2 + I6Q2 + (l9pi2 
+6pi3 + 12^23 ~ 27Qi + 2Q2)Q3 — I6Q3) + m^(l7p23Qi ~ 23p23Qi + 4Qi 
+ 4 Pi 3 (Qi ~ 5Q2) — 8P23Q2 — 4P23Q1Q2 + I2P23Q2 ~ 8Q1Q2 — 4Q2 + QpuQs 
+ (47'23 ~ 25P23Q1 + 12Qi + 4(^23 + Qi)Q2 ~ 4Q2)Q3 + 2( — 5^23 + 6Ql + 2Q2) 
xQg + 4Q3 +pi2(4pi3Qi + IIP23Q1 ~ 4Q3 — I3P13Q2 — IOP23Q2 + 8Q1Q2 
+8Ql + (I6pi3 + 13P23 - 13Qi - 8Q2)Q3 - IIQ3) +Pi3(20p23Qi - 8Ql 
—29P23Q2 + 23Q1Q2 + 24Q2 + QP23Q3 ~ I2Q1Q3 + IIQ2Q3 — 8Q3)) + (2m® 

+3 (pi3 + P23 — Ql — Q2) (P23Q1 ~ P13Q2) + 2pi2 (p23 ~ Q3) + (4pi2 + Pl3 

+2^23 + 6Q1 — IOQ2 + 9Q3) +Pl2(4p23 +^13(4^23 ~ 2Q2 — Q3) + 3 (Qi — Q2) 
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xQs ~ P 23 ( 2 Qi + 4Q2 + 5Q3)) + m'^{pi 2 — Pi 3 + 9P23Q1 ~ — QP23Q2 

+3Q1Q2 + 3Q2 + 5 ^ 23^3 ~ QQ1Q3 ~ 3Q3 + ^12(2^13 + 8 p 23 — 2 (( 5 l + 5Q2) 

+Q3) + Pl 3 {p 23 + 8 Qi — I3Q2 + QQs) ) ) + 4 ( — (pi3 + P23 — Ql — Q2) { — P23Q1 
+P13Q2) {pi 3 + P23 — Ql — Q2 — Q3) + P12Q3 {pi 3 + Ql + Q3) + P12 (Pi3 { — Q2 
+Q3) + (P 23 ~ Ql ~ Q2 ~ Q3) (P23Q1 + (Ql ~ Q2) Q3) + Pl 3 (Q2 (Ql + Q2 ~ Q3) 
+P 23 (Qi ~ Q2 + Q3))))) 5 ( 5 . 23 ) 

= 2 M(l 0 m® + 4 (pi 2 +P13 +P23 - Qi - Q2) (^23(^13 +P23 - Qi) +P13Q2) + 

X ( 4 M^ + 7 pi 2 + 12 pi 3 + 23 p 23 ~ 15 Qi + 7Q2 — 2IQ3) + 2 ( 2 p ^2 + Pi 2 ( 2 pi 3 + P23 
— 2 Qi) — 2^23(^13 + P23 ~ Ql) ~ (pi 3 ~ P 23 )Q 2 )Q 3 ~ 2P12Q3 + m^( 4 p ^3 + 20 pi 3 

xp 23 + 17^23 ~ SPlsQl ~ 23P23Q1 + 4Q3 + I4P13Q2 — 3P23Q2 — 4Q1Q2 — 7Q2 
+Pl2(4pi3 + 11 P 23 - 4 Qi + IIQ2 - 4Q3) + M^( - 2pi2 + Pi3 + 5 p 23 “ Ql + 2Q2 
~ 2 Q 3 ) — I3P13Q3 ~ I9P23Q3 + I2Q1Q3 + 5Q3) + {p 23 {pi 3 + P23 — Ql) 

+P13Q2+Pi2( - 2 p 23 + Q3))) , ( 5 . 24 ) 

aifu = - 2 M(l 5 m® + 4 (pi 3 +p 23 Ql+Pl 3 (p 23 -Q 2 ))(pi 2 +Pl 3 +P 23 -Ql “ Q2) 
~ 2 (Pi 2 + 2 Pi 3 + P23Q1 + Pl 3 ( 2 p 23 ~ Ql ~ 2Q2) + Pl 2 (pi 3 ~ Ql + Q 2 ))Q 3 
-2P12Q3 + [ 3 rri^ + pl^ + P23 {^pu + Qi) + P13{p 23 - Q2) + { 3 pi 2 + 4pi3 

+ 3 p 23 + Ql ~ Q2 + Qs)) + ? 71 ^(lOpi 2 + 31 pi 3 + 23 p 23 ~ 3 Ql — 2(8Q2 + Q3)) 
+m^(20p^3 + 8^23 + P23Q1 ~ 4Q3 — I2P23Q2 + Q1Q2 + 4Q2 — 4P23Q3 ~ 5Q1Q3 
+Q2Q3 ~ 3Q3 +Pi 2 (l 3 pi 3 + 10 p 23 + 3 Qi — 8Q2 + Q3) +P13 ( 29^23 ~ 9 Qi 
- 8 ( 3 Q 2 + Q 3 )))), ( 5 . 25 ) 

5 ' 2 sm = 2 M( 8 m® + m‘^(2M^ + 19 pi 2 + 6^13 + 12^23 - 6Q1 - I6Q3) + m^( 9 pi 2 + 6^13^23 
'^^P 23 + P13Q1 ~ 6P23Q1 + 3P13Q2 — Q1Q2 ~ 3Q2 + Pl2(l6pi3 + 13 p 23 ~ 9 Qi 
—7Q2 — IIQ3) + M‘^{ 2 pi 2 — Pi 3 + 2 p 23 + Q2 — Q3) ~ 8P13Q3 — IOP23Q3 + 7 Qi 
XQ3 + Q2Q3 + 4Q3) + Pl2(^^(pi3 + P23 ~ Q2) + 2(2pi2Pl3 + 2p^3 + 2 pi 3 P 23 
—2pi3Qi — P23Q1 — P12Q2 — 3P13Q2 — 2P23Q2 + Q1Q2 + Q2 + (2pi2 ~ Ql + Q2) 
XQ3))), ( 5 . 26 ) 

fsM^ ~ 4 (m®( 7 (Qi + Q2) + 4Q3) + m^( — (Qi + Q2) (l 2 pi 2 — llpi 3 ~ IIP23 + 4 (Qi 
+Q2)) ~ 2 ( 7 pi 2 — 3 pi 3 — 3^23 + 2 (Ql + Q2)) Q3) + 2pi2 { — 2 (pi3 + ^23) (Pl2 
+P13 + P23 - Ql - Q2) (Ql + Q2) + (pi2 - Pl 3 - P23) {Pl2 + Pl 3 + P23 - 2 (Ql 
+Q2))Q3 ~ 2P12Q3) + m^( — P 32 ( 10 (Qi + Q2) + 7Q3) + 2(pi3 + P23) (2(pi 3 
+P23 ~ Ql ~ Q2) (Ql + Q2) + (pi 3 + P23 ~ 2 (Qi + Q 2 ))Q 3 ) + Pl2( ~ 13 pi 3 (Qi 
+Q2) ~ 13 p 23 (Ql + Q2) — 12 (pi3 + P23) Q3 + 8 (Ql + Q2 + Q3) ) ) + 

X ( - 4 m® + m"^ (8pi2 + 3 ( - 2pi3 - 2 p 23 + Qi + Q2)) + P12 ((7*13 + P23) (27*12 
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+ 2 pi 3 + 2 p 23 ~ 3 (Qi + Q2)) + 3P12Q3) + ^^(6pi2 — {pi 3 + P23) ( 2 pi 3 + 2 p 23 
— 3 (( 5 i + Q2)) + ^12(8^13 + 8 p 23 — Q{Qi + Q2 + Q3))))) ) (3-27) 


/2442 

‘SM 


^ 2442 

92 SM 


AM ( - 7m® + P12 (pi3 + P 23 ) + 4 (pi2 + P13 + P23 - Qi - Q 2 )) + 2pi2 {pu 

-Pl3 - P23)Q3 + m'^{- + I2pi2 - llpi3 - llP23 + A{Qi + Q 2 ) + 2 Q 3 ) 

+m^(l0p^2 + ^^(2pi2 ~ P13 ~ P 23 ) + 2(pi3 + P 23 ) ( ~ 2pi3 — 2p23 + 2(^Qi + Q 2 ) 


+ Pl2(l3pi3 + 13p23 ~ 8 (Qi + Q 2 + Q 3 )))) ; 

(5.28) 

- 4m® + 2m'^(7pi2 - 3pi3 - 3p23 + Qi + Q 2 ) +Pi 2 {- Af^Pi2 - 2(pi2 

-Pl3 

(pi2 + Pi3 + P 23 — Qi — Q 2 ) + 4P12Q3) + m? {? 2 M‘^pi 2 + 7^32 ~ 2 (pi 3 

+ P23) 

+ P 23 -Qi- Q 2 ) + 4pi2(3pi3 

+ 3p23 ~ 2 (Qi + Q 2 + Q 3 )))) , 

(5.29) 

_ 1331 7I33II 

i/lSM “ ‘SM |piop2 ’ 

^ 2442 72442 

yiSAd — '-SAd 5 

(5.30) 

r2332 _ ^14411 

J SAd J SAd ^ 

72332 ^ 14411 

'SAI — 9lSAd ’ 

(5.31) 

^ 2332 714411 

9lSAd — '-SAd |pio-p2 ’ 

„ 2332 _ 14411 
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SM y^SM ’ 

(5.32) 

3443 _ 712211 

y^SM SM \p^^p2_^rni^m ’ 

^ 3443 _ ^ 12211 

y^^SM y‘^SM 



Appendix C 

For decays with five or more particles in the final state, not all of the scalar products of 
their momenta can be written as linear combinations of the invariant variables. In these 
cases, symmetry considerations allow to derive their expressions [24, 52], Although the 
procedure is quite transparent, results are not always compact. In our case this is so for 
a particular product, p^ ■ p^, whose explicit expression we quote in this Appendix for the 
reader’s convenience. 


P3-PA 


(m^ + - U2) (53(53 + t 2 - ml) (mi + U2 - 2 m^) - 253(^2 - ^3) (^i + M2 

-2m^) - 2^3^2(51 + Ml + M2 - 2m^ - - ml) + (S3 + ^2 - mf) (^2 - ^3) (5i 

+Mi + M2 - 2m^ - - ml) + 453^2(Af^ - M3) - (53 + ^2 - ml)‘^[M'^ - M3)) 

+ (52 - 53 - m^) (4M^S3t2 - 2 M^(s3 + ^2 - 'ml) (t 2 - ^3) - 53 (mi + M2 - 2m^)^ 
+ (^2 - h) (wi + M2 - 2m^) (si + Ml + M2 - 2m^ - - mf) + (53 + ^2 - mt-i) 

X (mi + M2 — 2m^) (M^ — M3) — 2^2(51 + Ml + M2 — 2m^ — — ml) (M^ — M3)) 

+ (mi -t 2 - m^) { 2 M‘^s^{ml - 53 - ^2) + 4M^S3(t2 - ^3) + 53(mi + M2 - 2m^) 

X (si + Ml + M2 — 2 m? — — m^) — (^2 — ^3) (51 + Mi + M2 — 2m^ — 

-m^)^ - 2 s3(mi + M2 - 2 m^) (M^ - M3) + (53 + ^2 - m?^ (si + mi + M2 - 2 m? 
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- ml) (M^ - M3)) 4 ( 4 M^S 3 t 2 - M^(s 3 + t2 - - ssiui + U2 

—2m?Y + {^3 + t2 — ml) (mi + M 2 — 2m^) (si + Mi + M 2 — 2m^ — — ml) 


■^2(51 + Ml + M2 — 2 m^ — — Mii)^) 


-1 


( 5 . 33 ) 


Appendix D 

Here we give the expressions for '^rl^ '^lrrl '^llrr contributions introduced in 
Section 3 in terms of the reduced form factors /, gi, g2 and I, which are dehned in the 
same way as those of the SM contribution (Tsm)- 


(i) For the contribution, in terms of the reduced form factors (we use the generic 
expression x = f, I, gi, 5^2), we have 


5,11 

^RL ~ 


^11 

■^SM 

^+22 

M2 

■^SM 

^S,33 

^33 

'^SM 

4 ’ 

"'"ML ~ 

4 ’ 

•^RL ~ 

4 

SM 

^RL — 

^44 

■^SM 

4 ’ 

S,2442 
^RL — 

™2442 

■‘^SM 

4 ’ 



( 5 . 34 ) 


^5,1221 

Jrl 


( 3 m^(pi2 + P13) + 3 (pi2 + ^13)^23) +ml{- 3 m^{Q2 + Q3) - 3 p 23 

X [Q2 + Qs)) + 4 m^(pi 2 (Qi + 2Q2) ~ Qi{ ~ P23 + Qi + Q2 + Qs) + Pi 3 
X (Qi + 2Q3)) + 4 m^(p^2 {Qi + 2Q2 ~ Q3) + Qi{~ P23 + Qi + Q2 + Q3) 
X ( — 2 p 23 + 2Qi + Q2 + Qs) + plsiQi ~ Q2 + ‘^Qs) ~ Pisi^Ql + Q1Q2 
—Q\ + P23 ( ~ 3 Qi + Q2 ~ 3Q3) + 6Q1Q3 + 3Q2Q3 + 2Q3) + pi2 [3P23Q1 
— 3 Ql + 3P23Q2 — QQ1Q2 — 2Q2 ~ {p 23 + Ql + ^Q 2 )Q 3 + Q3 + P 13 ( 2 Qi 
+3 [Q2 + Q3) ) ) ) + 4 ( — P13Q2 + P23Q1 {p 23 — 2(^1 — Q2 — Qs) {p 23 — Ql 
—Q2 — Q3) ~ P12Q3 + Pl 3 (p 23 {Ql ~ 2Q2 + Q3) + Q2 (3Q1 + 2Q2 + Q3) ) 
^ _|_ Q2 — Q3) {Ql + Q2 + Q3) ~ P23 

■ Qs) (3Q1 + 2Q2 + 


+^13(^23(2^1 ~ Q2 + Qs) ~ Q 2 ( 2 Qi + Q2 

Qs )) + P12 (P13Q2 + P23 {Ql + Q2 — 2Q3) 
2 ^ ^ \ . 2 ^ ^ 


X {Qi — Q2 + Qs) (3Q1 + 2Q2 + Q3) ) + P12 (P13Q2 + P23 {Qi + Q2 — 2 i 
+Q3 ( 3 ( 5 i + Q2 + 2Q3)) + P12 (p 23 (2Q1 + Q2 — Qs) + pIsQs ~ Qs (2Q 
—Q2 + Qs) {Ql + Q2 + Qs) ~ P 23 {Qi + Q2 ~ Qs) (3Q1 + Q2 + 2Q3) 4 ^ 

X (Q2(Qi + Q2) + Q1Q3 + Qs + P 23 ( 2 Qi + Q2 + Qs)))) + M"^ {m^{ a{pi2 

+P13) + 2Qi) + Pi3(4p23 ~ 5Q2) + P12 (4p23 ~ 5Q3) + ml { — 4m^ + P12P23 
+P13P23 + 4p23Ql - 4P13Q2 - P23Q2 + 'm‘^{pi2 + Pl 3 - 4p23 + 4Qi - Q2 
“Qs) ~ (4pi2 + P23) Qs) + P23Q1 (4p23 ~ 3 {2Q1 + Q2 + Q3) ) + (6pi2 

+ ^13(8^23 Ql — 6Q2) + Pi2(8pi3 + 8p23 Ql — 6Q3) + 3Qi { 2 p 
— 2 Qi — Q2 — Q3) ) + P12 (4^23 + P23 {Ql ~ 2Q2 ~ 8Q3) + 3Q3 (2Q1 — Q- 
+Q3) + Pl3 (4p23 ~ 3 {Q2 + Q3) ) ) + Pl3 (4^23 + 3Q2 (2Q1 + Q2 — Qs) 
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75,1221 

^RL 


+P 23 (Qi ~ 2 ( 4 Q 2 + Qs)))) + ^l( ~ ^p\^Q2 — + ^^( 3 pi 2 Ql 

+ 3 pi 3 Qi + 6P23Q1 ~ ^Qi ~ QP13Q2 — 8P23Q2 — Q1Q2 + QQ2 ~ (6pi2 

+ 8^23 Ql — 8 Q 2 )Q 3 + 6Q3) + 2 m'^(( 5 l — 2 (Q 2 + Qs)) + P 23 ( ~ ^Ql 
+ 4 p 23 (Ql — Q2 — Qs) ~ Qi{Q2 + Qs) + 4 (Q 2 + Q2Q3 + Q3)) + ^13(^2 
X (6Q1 + 5Q2 + 3(53) + P23 (3Q1 — 2 (4Q2 + Qs ))) + P12 ( 3 pi 3 [Q2 + Q3) 


+Q3 {QQi + 3(^2 + 5Q3) + P23 (3Q1 — 2 (Q2 + 4Q3)))) 


( 5 . 35 ) 


M 

y 


4771 “^ (pi2 + Pl 3 + P23 - Ql) + + 2P23Q1 + {Pl2 + Pl 3 


+ 6 p 23 + 2 Qi) + Pi 3 (p 23 - 2Q2) + P12 {p 23 “ 2Q3) ) + "^1 ( “ “ 3^13 

X (P 23 + 2Q2) + 2 p 23 ( ~ 2 p 23 + 3 Ql + Q2 + QQ) — 3 pi 2 (^23 + 2Q3) + 

X ( — 3pi2 — 3pi3 + 2( — 3p23 + 3(^1 + Q 2 + Q 3 ))) ~ 2m^(2p^2 + 2pi3 
+Pl 3 ( 6 p 23 - QQl +Q2- 4Q3) + Pl 2 ( 4 pi 3 + Qp 23 “ 6Q1 - 4 Q 2 + Q3) 

+2 ( 2^23 + ‘^Qi + Q2 + Q2Q3 + Q3 + 2 Qi (<52 + Q3) ~ P23 (4Q1 + Q2 
+Q 3 ))) + 2( — 2pl^{p23 + 2 ( 52 ) — 2 ^ 12(^23 + 2 Q 3 ) — P23{‘2p 23 + ( 2 Q 1 
+Q2) + 4Q1Q3 + Q3 — 2 p 23 ( 3 Ql + Q2 + Q3)) + Pl 3 ( ~ 4^23 + Q 2 ( 4 Qi 


+3Q2 + Q3) + P23 {QQi — 2Q2 + 3(53) ) + Pi2 ( — 4^23 + P23 {QQi + 3(^2 


“ 2 ( 53 ) “ 4^13(^23 + Q2 + Q3) + ^3(4(51 + Q2 + 3(53))) ) 


( 5 . 36 ) 


5,1221 

9irl 


— 4pi3 + mi(3m^ + 3^23) + - 8pi2 + 4(5i) + 2p23Qi{2p23 - 2Qi 

-Q2 - 2Qs) - 2pi2(2(pi3 + P23) + Q3) + m\{Am‘^ - ?>pi 2 Pi 3 + 3pi3 
+2pi2P23 + 8pi3P23 + 4^23 “ P 23 Q 1 “ ^Pl3Q2 “ 4^23^2 + ^77^(6^13 + 8^23 
—Ql — QQ2 — 4(^3) “ (2pi2 + P13 + 2^23)Q3) + (p^3 + ml {m? + P23) 

-P 12 (P 13 + 2 ,p 23 ) + P 23 ( 5 i + m^( - 6 pi 2 “ 2 pi 3 + (^i) + ^>13(2^23 “ Q2 
+Q3) ) + Pl 3 ( 4^23 “ 8P23Q2 + ^QiQ 2 + 2 {Q2 — QQ) {2Q2 + QQ) ) + P13 

X (8p23 “ 2{2Qi + 4(^2 + Q3)) + 2pi2( — 2^23 “ 2^13(^23 — Ql + Q2) 
+^3(^1 “ Q2 + Q3) + P23{5Qi + 2((52 + Q3))) + 2m^( — 4p^2 + 2pi3 
“Qi( “ 4 p 23 + 2 Qi + (^3) + ^13(2^23 + 3 Qi — 2Q2 + 2Q3) + pi 2 ( — 6 pi 3 


“ 6 p 23 + 4 ( 2 Qi -h Q2 -h Q3))) , 


( 5 . 37 ) 


^ 5,1331 

JRL 


1 

4 


M^( 4 pi 3 P 23 - 3 p 23 Qi - 3P13Q2 - 3P12Q3) + + lOpuQi + 23 p 23 Qi 


— 4 QI + 23P13Q2 + IOP23Q2 “ 48 ( 5 i (52 “ 4 QI + 2 pi2(Qi + (52) + 2 M‘^ (Qpi2 
+ llpi 3 + llp 23 “ Q{Qi + Q2)) — 3IP12Q3 + 2 ( 5 pi 3 + 5^23 “ 11 (Ql + Q 2 ))Q 3 
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-'23 


—8Q3) + 772 .® (l6Af^ + 13 (Qi + Q2 + Qs)) + ~ 2 Qi ( 2 pi 2 ~ 7 Pl 3 P 23 ~ 5 p 2 

+P 12 (2pi3 + P 23 ~ 4Qi) — 2P13Q1 + 2(5i) — 2 (2p^2 ~ P13 (5pi3 + ^^23) + 20 (pis 
+P23 )Qi ~ + Pl2(pi3 + 2p23 + 4Qi))Q 2 + 4(2pi2 + P 23 + 4Qi)Q2 ~ 4Q2 

+ +f^( 4 pi 3 + 4 p 23 ~ 3 (Qi + Q2 ~ Q3)) ~ 2 ( 3 p ^2 ~ 2 (pi 3 P 23 ~ 5P23Q1 ~ 5P13Q2 
+ I3Q1Q2) + ^12(7^13 + 7p23 — ^{Ql + Q 2 )))Q 3 + 4 ( 5 pi 2 + Ql + Q2)Q\ 

+ 2 M^( 3 P 33 + 12 pi 3 P 23 + 3 P 23 “ 5 ^ 13^1 “ 8 ^ 23^1 + ^QI “ 8P13Q2 “ 5 ^ 23^2 

+2Q1Q2 + 3Q2 + ^12(5^13 + 5 p 23 ~ 4 (Qi + Q2) — 5Q3) — 2 (pi 3 + P23 + Qi 
+^2)^3 ~ 2(53)) + 2 M'^ ( 2 p 1 ^(^P 23 — Q2) +P 23 Qi{ ~ 2^23 + 3 Ql + Q2) — Pl2 
XQ3 + QQ1Q2Q3 + ^13(2^23 + ^2(^1 + 3Q2) ~ 2 p 23 ( 2 (Qi + Q2) + Q3)) + P12 
X (^13(4^23 ~ 3Q2 ~ 2(53) — P23(3 Qi + 2Q3) + Q 3 (Qi + Q2 + 4Q3))) + 4 (p ^3 
X<52(^23 — Ql + Q2 — Q3) + P12Q3 + Pl2( ~ P23Q1 ~ P13Q2 + (pi 3 + P23 ~ 2 [Qi 
+Q2) ) Q3) + P12 ( — P23 {pi 3 + P23 — 2 Ql) Ql — Pl 3 [pi 3 + P23) Q2 + 2P13Q2 
~{{p 23 — Ql)Ql + {p 23 ~ ^Qi)Q2 — Q2 + Pl 3 {p 23 + Ql + Q 2 ))Q 3 + (Pl 3 + P23) 

xQI — Ql) + ^13(^23^1 + P 23 {Q\ ~ 6Q1Q2 + Q2 “ (Ql + Q 2 )Q 3 ) + Q2(Qi 

+Q1Q2 ~ Q2 + 6Q1Q3 + Q3)) + Qi(p 23 (Qi ~ Q2 ~ Qs) ~ 4Q2Q3 (Qi + Q2 


+Q3) + P 23 ( ~ Ql + Q1Q2 + Q2 + 6Q2Q3 + Q3))) 


( 5 . 38 ) 


75,1331 

'-RL 


— - 13?71® + M^P23 + 771^ (4M^ - 2pi2 - 10pi3 - 23p23 + 4Qi + 24Q2 

+2Q3) + 2M^( — 3P23Q1 +Pl3(2p23 + Q2) ~ {p 23 + Q2)Q3 + P12 (P23 
+Q3)) + rn^ + 2M^(3pi2 + 4pi3 + 2 p 23 — 3Qi + Q2 — 3Q3) + 2(2p^2 
~^P23 + 2Qi + IOP23Q2 ~ 3Q1Q2 — 5Q2 + ^12(2^13 + P23 ~ 4Qi + 2Q2 
“Qs) + 4P23Q3 + 3Q1Q3 — 8Q2Q3 + 2Q3 — ^13(7^23 + 2Qi — IOQ2 
+3Q3))) + 4 (P13Q2 + P23Q1 ( “ P23 + Ql + Q2) + (p 23 + P23 (Ql “ 2Q2) 
+Q2 (2Q1 + Q2) ) Qs + Q2Q3 + P12 {p 23 + Qs) + P12 (p 23 “ P23 (2Q1 + Q2 
“Qs) + Pl3 (P23 + Q2 + Qs) “ Qs (Ql + 2Q2 + Qs) ) “ Pl3 (P23 + ^23 (Ql 


“3Q2) + Q 2 (Qi + 2Q2 + 3Q3))) , 


( 5 . 39 ) 


5,1331 
92 RL 


— - 13771 ® + -P12) + 771 ^(llpi 2 - 10 (pi 3 +P23 - 2 (Qi + Q2)) 

+8Q3) + 2771 ^ (3p^2 + 6P23Q1 “ 3Q3 + 3P23Q2 — IOQ1Q2 — 3Q2 — ^13(2^23 
“ 3 Qi — 6Q2) + ^12(7^13 + 7p23 — 3 (Qi + Q2) — IOQ3) — 4 (Qi + Q 2 )Q 3 ) 
+4 ( “ Pi2 + (pi3 “ Ql) {p23 “ Q2) (Ql + Q2) + P12 ( “ Pl3 “ P23 + Ql + Q2) 
“ (P23Q1 + (pi 3 “ 2 Qi) Q2) Qs + P12 ( “ Q1Q2 + Pl 3 {p23 “ Qs) + Qs (Ql 
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(5.40) 


+Q2 + Qs ~ P23))) + — 2 (p ^2 + P23Q1 ~ ^13(2^23 ~ Q2)) + 

X (l0pi2 + 8pi3 + 8 p 23 ~ ‘^{Qi + Q2 + 2Q3))) , 


^ 5,1441 

JRL 


- - 3 M'^pi 2 P 23 + - Qi + 14Q2 - 7Q3) - m‘^{Qpi 2 Qi + PisQi 

+5^23^1 ~ ^Ql ~ P12Q2 — 25P13Q2 ~ 22P23Q2 + I4Q1Q2 + I8Q2 + (l9pi2 

+9pi3 + llp23 ~ 19(5l + 4 ( 52 ) Q 3 ~ 12(^3 + M^(6pi2 — Pl3 — 2^23 + 3Qi 

-9Q2 + 8Q3)) - m^(3MVi2 + 3 pi 3 P 23 Qi + GplsQi - 4pi3Q? - IIP23Q? 
+ 4 Qi — IIP13Q2 ~ 22P13P23Q2 ~ 8 ^ 23^2 + 19 pi 3 ( 5 lQ 2 + ‘2P23Q1Q2 + 22 
XP13Q2 4 " 12 ^ 23^2 ~ 8Q1Q2 — 4(^2 + ( 2 pi 3 + 4^23 ~ 15 ^ 23^1 + 12 Q^ + 2 

XP23Q2 + 4 :QiQ2 — 4Q2 + ^13(8^23 ~ I6Q1 + 7Q2))Q3 + 2( — 5pi3 — 4^23 

+6Q1 + 2Q2) Ql + 4Q3 + P32 (4Q1 + 6Q2 + 4(53) + P12 (2 (2 (pi3 + 3p23 
—2Qi)Qi + (pi3 + P23 — 4:QijQ2 + Q2) + (l3pi3 + 10^23 ~ 17Qi — I6Q2) 
XQ3 — I3Q3) + iV 7 ^( 5 p 32 4 “ Pi 3 4 “ 5P23Q1 ~ 3Q3 — 5P23Q2 4 " 3Q1Q2 4 “ 3Q2 
~Pl3{P23 ~ 2Ql + IIQ2 ~ 7Q3) + 4^23^3 ~ 6Q1Q3 — 3Ql +Pl2(4pi3 + 10 
xp23 — Ql — 7Q2 + Q3))) 4- — P23Q1 4- P13Q2) (2pi3 + 2p23 — 3((5l 

+Q2) ) 4- P12 ( — 4p23 4- 3Q3) + P12 ( — 4^23 4- 3 ( — Qi+Q 2)Q3~ Pl3 (2p23 
—Q2 — 2Q3) + P23 (Qi 4- 2Q2 4- SQs))) 4- 2(2^32^3 4- ( — P23Q1 4- P13Q2) 
X ((pi3 4- P23 — Ql — Q2) {p23 — 2(^Qi + Q2)) — [pi3 + P23 — 2(^Qi + Q2)) 
XQ3) “^12(2^23(^1 4- Q2) 4- 2pi3(Q2 — Qs) — 3^23^3 4- Q3(4 Qi + Q2 
4-2(53)) + P12 ( — 2p\^Q2 + P23 ( “ 3(5 i — 2Q2 + Qs) 4-2 ((5i — Q2) Qs {Ql 
+Q2 4- Q3) 4- P23(4(5 i 4- 7Q1Q2 + 2 Q\ — QxQj, -|- 3Q2Q3 ~ QT) ~ Pi 3 {Q 2 


+P23{2 ‘Qi + 3Q2 — Qs) -|- 2Q1Q3 + (^s))) ) 


(5.41) 


, 5,1441 

''RL 


— + 771^{- 3M^ + 6pi2 -h Pi3 5p23 - 5(5i IIQ2 - 13(53) -4 

X ( - P23ipi2 -4 Ql) + P13Q2 -4 P12Q3) + "i^(4pi2 -4 3pi3P23 + 6^23 “ 4pi3(5i 

— IIP23Q1 4- 4(5^ -4 14^13(52 4- P2SQ2 — 2QxQ2 — 3Q2 + ^12(4^13 -4 12^23 
—8Qi + 3Q2 — QQs) ~ IIP13Q3 “ "^P2sQs 4- 10(5i(53 4- 2Q2Q3 + 7Ql — 

X (pi2 -4 3^23 + Ql — 2Q2 + 2Q3)) + 2(2^33(52 4- 2^32(^23 4- Qs) 4- P23{p 23 
—8p2sQi + 2QI — 2P23Q2 4- 3Q1Q2 + Q2 + QiQs) 4- ^12(3^23 “ Q3{2‘Qi 
-4(53) 4- P2S ( “ 4(5i — 3Q2 + Qs) + 2,pi3 {p23 + Q2 + Qs) ) 4- Pl3 (P23 "4 P2S 


X ( — 2Qi + Q2) — Q2{2 ‘{Qi + Q2) 4- Qs))) ; 


(5.42) 
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S,1441 

9irl 


14 m® + M^(pi 2 (pi 3 + 2^23) - P23Q1 + P13Q2) - ^"^(^12 + 25^13 
— 3 (( 5 i + 6Q2) ~ 2Q3) + m^(6p^2 ~ llPi 3 ~ 22 ^ 13^23 


8 pL 


M 

T 

+ 22^23 

+ 5 pi 3 Qi + P 23 Q 1 + 2 Q^ + 22 P 13 Q 2 + I 2 P 23 Q 2 ~ 5 Q 1 Q 2 ~ 4 Q 2 + P 12 
X ( 2 pi 3 + 2^23 ~ llQi + 2Q2 ~ 9Q3) + M^( 4 pi 2 — Qi + Q2 ~ Q3) 

+ ( 4 pi 3 + 2 p 23 + 3 Qi — 3 Q 2 )Q 3 + QI) + 2 ( 2 p^ 2 (pi 3 +P23) ~ P 13 {p 23 

— 2Q2) — P23Q1 {p 23 — Ql — Q2 — Q3) — Pl 3 {p23 ~ 3 ^ 23^2 + Q2 (2 [Qi 
+Q2) + Qs)) +^12(2^13 + 2^23 +^13(3^23 ~ 2 Qi + Q2 — Q3) + Q2Q3 


■P 23 ( 4 Ql + 2Q2 + 3Q3))) 


( 5 . 43 ) 


S,1441 

92 RL 


— — - 7 m® + M‘^pi2{pi2 + 2 p 23 + Q2) + - 19 pi 2 - 9pi3 - llp 23 

+6Q1 — 2Q2 + 12(53) + 2 ( 2^32 + Pl2 {p 23 ~ Q2) {pi 3 + P23 + Q2) — {pi 3 
+P23 — Ql — Q2) {P13Q2 — P23Q1) + Pl2(2pi3 + 3 p 23 ~ 2 Qi — Q2 — 2 Q‘^ 
-Pl 2 {pi 3 + P23 -Ql + 2Q2)Q3) +m^[- Apl^ - 2^33 - %Pi 3 P 23 “ 4^23 
+^P 13 Qi + 8P23Q1 ~ 2 Q\ — 3^13(52 + Q1Q2 + Q2 + (l0pi3 + 8^23 ~ 5 ( 5 l 
—3Q2) Qs — 4(^3 + ( — 2pi2 — 2pi3 + P23 — Q2 + Q3) + P12 ( — 13 pi 3 


— 10^23 + llQl + 7(^2 + 13 ( 53 ))) 5 


( 5 . 44 ) 


S,1221 _ 5,12211 

92RL — 91rl \p2-H-p3 ’ 


5,1331 _ 75,13311 
91rl RL Ipi-f^p2 ’ 


( 5 . 45 ) 


rS,2332 _ ^5,14411 
" RL " RL ^ 

5,2332 _ ,5,14411 
91rl RL |pj4_^p2 ’ 


75,2332 _ 5,14411 

RL 91rl \p^j^p2 ’ 

5,2332 _ 5,14411 

92rl — 92rl ’ 


( 5 . 46 ) 


r5,3443 
J RL 


fSp221 

J RL 




^5,3443 _ 5,12211 

RL 91RL \p^j^p2^rni^m ’ 


5,3443 

9irl 


75,1221 

'-RL 


pi<^P2,mi^m ’ 


5,3443 

92rl 


5,1221 
92 RL 


pi^P2,mi^m * 


(5.47) 


{ii) Regarding the expressions for the TR contribntion, the reduced form factors 
shown below must be multiplied by the common factor (M^ + m^ + 2 (m^ +P12 +P13 +P23 
—Ql — Q2 — Q3)) (iLii = m if there are identical particles): 


Irl^ = 12 [M^(m^ + P 23 )(Pi 2 + P13 ~ Ql) ~ Ln^Qi + P12P23Q2 — P23Q1Q2 + P13Q2 
+P13P23Q3 ~ P23Q1Q3 — P12Q2Q3 — P13Q2Q3 + 2Q1Q2Q3 + P12Q3 + Ln? 

X ( —Qi(P23 + Q2 + Q3) + P12{2Q2 + Q3) + Pl3(Q2 + 2 Q^))] , ( 5 . 48 ) 


= 12 M [(m^ + M^){m^+P 23 ) - 2Q2Q3] , 


( 5 . 49 ) 
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= 12 M [P23Q1 - P12P23 + m^iQi - 2 pi 2 - P13) - P13Q2 + P13Q3] , ( 5 . 50 ) 

^ = 12 \2p12P13Q1 + P12P23Q1 + P13P23Q1 + P12P13Q2 - m^Qi - PI3Q2 + P12P23Q2 
“^ 12^3 + P12P13Q3 + P13P23Q3 ~ m\{nn? + P 23 )(Qi + Q2 + Q3) + ^^{—P 23 Qi 
+P12(Qi + 2(^2 + Q3) + P 13 (Qi + Q2 + 2(53))] ) ( 5 . 51 ) 

~ 12M + 2(m^ + pi2 + P13 + P23 — Qi — Q2 — Q3)] 5 ( 5 . 52 ) 

= 12 [M^(Pi 2 +Pl 3 )P 23 - 2 m^Qi + P12P23Q1 + P13P23Q1 - ‘^P 23 QI +P12P13Q2 
—P13Q2 ~ f^\P 23 Q 2 + 2P12P23Q2 + 2P13Q1Q2 ~ P23Q1Q2 + P13Q2 ~ (Pl2 + P23 
X(m^ — 2^13 + Qi) + P13Q2 + ^12(^2 ~ 2 ( 5 i — Pi 3 ))Q 3 + P12Q3 + 

X (pi2 + P13) + P13Q1 ~ 2 ^ 23^1 ~ 2 Q^ + 2P13Q2 — Q1Q2 + 4P13Q3 — Q1Q3 
—m ^((52 + Q3) + Pi 2 (Qi + 4Q2 + 2Q3))] , ( 5 . 53 ) 

= 12 ^ [ 2 m'^ - P13P23 + 2P23Q1 - rn^{pi2 + P13 - 2 (p 23 + Qi)) - 2P13Q2 

—Pi 2 {P 23 + 2Q3)] , ( 5 . 54 ) 

9 irl^^^ = 12 M [p 23 (’^? + Qi) - ^12(^13 + 2^23) + rn‘^{ml - Apu - 2pi3 + Qi) 

+Pi 3 (Pi 3 — Q2 + Q3)] ) ( 5 . 55 ) 

V,ll _ V,lli V ,22 _ ,V ,22 

92rl ~ 91rl Ip2 - H . p3 ) 91rl — ^RL 1 


1^,22 y,22| 

92rl — 9irl \p2^P3 5 

y,1221 

92 RL 

V,1221, 

~ 9irl P2'H-p3 7 

(5.56) 

rV,33 

JrL ~ JRL piOp2,mi->.m ; 

jV,33 
''RL — 

y,ii| 

91RL pi<->'P2,mi->.m 7 


V",33 _ 7^,11 1 

91rl ~ ^RL lpi«'P2,mi->-m) 

92rl — 

Kill 

92jil \pi-^P2,'nT'i^'m ? 

(5.57) 

rVM _ fV,22i 

JRL ~ JRL \pi<r^P 2 ,mi->-m y 

]V,44 
''RL ~ 

V,22\ 

91RL lpiOp 2 ,mi —7 


V,44 ,V,22 1 

91RL ~ ^RL piOp2,mi^.m ; 

V,44 

92rl — 

K22| 

92rl \pi-H-p 2 ,mi^m 7 

(5.58) 


~ —6(m^(2M^ — Qi — Q2 + Q3) + M‘^{P 23 Qi + P13Q2 + P12Q3) + 'm‘^{ 4 :QiQ 2 
—P23Q1 — 2 ,QI — P13Q2 — 2Q2 + 2pi2(Qi + Q2) + M‘^{2,pi2 + 2 pi 3 + 2 p 23 
+Ql + Q2 — Q3) + 5P12Q3 + 2 (Pi 3 + P23 ~ 2(Ql + Q 2 ))Q 3 ~ 4Q3) + 2 
X ((Q2 ~ Qi)(P 23 Qi ~ P13Q2) ~ P12Q5 ~ 2Q1Q2Q3 + P12(P23 Qi + P13Q2 
+ (Ql + Q2)Q3))) ; (5.59) 
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~ — 6 Af(m^ — M‘^P 23 — rri^(^M‘^ + 2 pi2 — P23 ~ 2Q\ + 2Q2 ~ ‘^Qs) 

+2,{p23Qi — Q 2 {P 13 ~ Q3) ~ Pl 2 {P 23 + Qs))) ) ( 5 . 60 ) 

~ 6M(m^ — pi2(M^ + 2^12) + m^(M^ + 5 pi 2 + 2 (pi 3 + P23 — 2 Q^))) , ( 5 . 61 ) 

= - 12 {M‘^pi 2 P 23 + P12P23Q1 - P23QI + P12P13Q2 + 2P12P23Q2 + P13Q1Q2 
—P23Q1Q2 + P13Q2 ~ + Q2) — PuiPu — Qi + Q 2 )Q 3 + m^{M'^pi2 

—P23Q1 ~ Qi ~ P13Q2 — P23Q2 + Q1Q2 + Q2 + P13Q3 ~ 2.Q1Q3 — QI 
+P 12 (Qi + 3 (Q 2 + Qs)))); ( 5 . 62 ) 

~ —12M(m^ + P23Q1 ~ P13Q2 — Pi2{P23 + Q3) + — 2Q2 + 2(^3 

-P12+P23)), ( 5 . 63 ) 

9 irl^^^ ~ — 12 M(m^ — ^12(^13 + 2^23) + P23Q1 ~ P13Q2 + rn^{Qi — Q2 + Q3 

- 3 pi 2 + Pi3 + P23)), ( 5 . 64 ) 


= - 12 M(pi 2 (pi 2 + Q2) - m^{ 3 pi 2 + P13 + Q2 - Q3)), 


( 5 . 65 ) 


yy^442 _ 24(^(^777,4 — Pi2(Pi3 + P 23 ) + + Pl3 + P23))(Ql + Q 2 ) 

-(2m^ -P12)P12Q3) , 

- 24 M(m^ - Pi2(pi3 + P 23 ) + m^(pi 3 + P 23 - 2^12)), 
- 24 Mpi 2 (pi 2 - 2m^), 

V,1331 _ 7y,133] 

91rl — ^RL 

rV,2332 _ fVpiAl 
JrL — JRL 

y,2332 _ ,y,1441| 

91rl — ^RL \ 


jV,2A42 

^RL 


y ,2442 

92 RL 


fV,34A3 _ fV,l221 
J RL ~ J RL 


> RL 
y ,3443 


9irl = I 


yi22i 


ll 

P1-H’P2 ) 

y,2442 

9irl 

,y,2442 1 

" ^RL \pi^P2 I 

\pi-^P2 1 

iV,2332 
^RL — 

y,1441| 

9irl \pi^P21 

Pl^P2 ) 

y,2332 

92 RL — 

y,i44i| 

92rl \pi^P2 

\pi<^P2 ) 

iV,3443 
^RL ~ 

y,i22i| 

9irl \pi-<r^P2 ) 

Pl^P2 ) 

y,3443 

92 RL — 

y,i22i| 

92rl \pi^P2 


( 5 . 66 ) 

( 5 . 67 ) 

( 5 . 68 ) 

( 5 . 69 ) 


( 5 . 70 ) 


( 5 . 71 ) 


(Hi) For the contribution all reduced spin-dependent form factors are zero. 

The results of all reduced spin-independent form factors, which are shown below, must be 
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multiplied by the common factor + m\ + 2(m^ +pi 2 +P 13 +P 23 ~Qi—Q 2 — Qz)) 

(mi = m if there are identical particles): 

fLRRL ~ —6144(2m^ + 3m^p23 + P 23 + + P 23 ) — 2 Q 2 Q 3 )) (5-72) 

fLRRL = -6144(2m‘^ - 2pi2Pi3 + 3mV23 + pIs + mKrri^ + P 23 )), (5.73) 

/Sf ^ = -12288(2m^ + P 23 + P 23 Q 1 + m‘^i^P23 + Qi) - P 13 Q 2 - P 12 Q 3 ), (5.74) 

rSV,33 _ /■syil I rSVM _ rSV,22i 

JlRRL ~ JLRRL\pi-<^P2,mi^m , JLRRL ~ JLRRL\pi*^P2,mi^m ■, 

fLRRL^ = -3072(M^(pi 2 + Pl3 + P 23 ) - 3m^ + W?{M‘^ + Pi 2 - 3pi3 - 3p23 
+ 2 Qi + 2 Q 2 — 6 Q 3 ) — 2((—P 12 + Ql + Q 2 )Q 3 + P23(~Qi + Q 3 ) 
+Pl3(2p23 — Q 2 + Q 3 ))) , 

~ 6144(2m'^ + 2 pi 3 P 23 + P 23 + P 23 Q 1 ~ 2 P 13 Q 2 — ^ 12 ( 2^23 + Q 2 ) 

+P 13 Q 3 + fn^{Qi — 2 Q 2 + 2 Q 3 — 2 pi 2 + 2 pi 3 + 3 ^ 23 ))) 

fLRRT'^ = 6144(3m^ - P 12 - 2pi2(Pl3 + P 23 ) + {Pl3 + P23f 
+4m^(pi3 +P 23 -P 12 )), 

rSV,2332 _ rSV,IU1\ ^Sn,3443 _ nSV,1221, 

JLRRL ~ JLRRL \pi^P2 i JLRRL ~ JLRRL \pi^P2,mi^m ■ 

{iv) As in the previous case, for the TI^rr contribution only the reduced spin- 
independent form factors are different from zero. The expressions of the reduced form 
factors {i = 1, 2, 3, 4) and given below have to be multiplied by the com¬ 

mon factor Mmi{M‘^ + m\ -|- 2 {w? + pi 2 + Pi 3 + P 23 — Qi — Q 2 — Q 3 )) ("^i = m ii there 
are identical particles): 

fLLRR = -6144(2m^ -h 3m^p23 + PI 3 + + P 23 ) - 2 Q 2 Q 3 ), (5.81) 

fLLRR = -6144(2m'^ - 2pi2Pi3 + 3m^p23 + P 23 + + P 23 )), (5.82) 

fLL’RR^ = 4096Mmi(—4m® — 4m'^pi2 — 2im?p\2 — 4m‘^pi3 ~ ‘2rri^Pi2Pi3 — 2m^pl^ 

-10mV23 - 6m^pi2P23 - p\2P23 - 6mVi3P23 - p\3P23 - 8mV23 - 2 pi2pL 
- 2 pi 3 pL - 2 p 23 - Sm^Qi - 6m^pi2Qi - bm^pisQi - 14 mV 3 Qi - 6P12P23Q1 
-QP13P23Q1 - 6P23Q1 + Qrn^Ql + QP23QI + 4 m '^(52 + 4mVi2Q2 + 8 m^pi 3 Q 2 
+6P12P13Q2 + 6P13Q2 + 6 m^p 23 Q 2 + 2P12P23Q2 + QP13P23Q2 + 2.P23Q2 
+QTn^QiQ2 - QP13Q1Q2 + QP23Q1Q2 - 2 m^Ql - QP13QI - P23QI + 4m^(53 


( 5 . 76 ) 

( 5 . 77 ) 

( 5 . 78 ) 

( 5 . 79 ) 

( 5 . 80 ) 
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^ 5 y ,1331 

Jllrr 


fSV,144.1 
JLLRR 


fSV,2442 
JLLRR 


+8mVi2Q3 + QPuQs + ^m'^PisQs + QP12P13Q3 + + 2P13P23Q3 

+2P23Q3 + Qm^QiQs — QpuQiQs + 6P23Q1Q3 ~ 2 m^Q 2 Q 3 ~ QP12Q2Q3 
-QP13Q2Q3 - 2 m^Ql - QpuQl - P23QI + + m^(p23 - 3Qi) 

+3 (Pi 3Q2 + P12Q3 ~ P23Q1)) + Ln\{rn^ + m^{p 23 ~ 3Qi) + 3(pi3Q2 

+P12Q3 ~ P23Q1)))) (5.83) 

rSV ,33 _ I rSV ,44 _ rSV, 22 \ re; S4'\ 

Jllrr ~ JLLRR\pi-<^P2,mi^m , Jllrr ~ JLLRR\pi-^P2,mi^m , 

-1024mM(llm® + 3 M^{pi 2 + P13 + P23) + m^{ 22 M'^ + 43pi2 + 29pi3 + 29p23 
—8Q1 — 8Q2 ~ 4Q3) + 2M^(6p^2 + 6pi3 + 6^23 ~ SpaaQi ~ 4^23^2 + Pi3(12p23 

—4Qi — 5Q2 — 4Q3) + 4pi2(3pi3 + 3p23 — Qi — Q2 — Qs) — 4P23Q3 ~ 3Q1Q3 
—3Q2Q3) + 2(4p]^2 + ^Pi 3 + 4p23 ~ 3P23Q1 ~ P23Q1 ~ 4P23Q2 ~ P23Q1Q2 + 2 Q^ 
xQa + 2Q1Q2 + 2 p^ 2 ( 6 pi 3 + 6 p 23 ~ 2 Qi — 2Q2 — SQ^) — QP23Q3 ~ I2P23Q1Q3 
+ 6 Q^Q 3 — 4P23Q2Q3 + 20 ( 5 iQ 2 Q 3 + 6Q2Q3 + 8Q1Q3 + 8Q2Q3 + ^13(6^23 
—4Qi — 3(Q2 + 2Q3)) + Pi2(12p^3 + 12p23 — 2Q1Q2 + Pl3(18p23 ~ 8Qi — 6Q2 
“QQa) ~ JQ1Q3 ~ JQ2Q3 ~ P23(6 Qi + 8Q2 + 9Q3)) + Pl3(6p23 ~ 3p23(Ql 
+Q2) ~ Qi(Q 2 + 4(^3) ~ ^2(^2 + I2Q3))) + rn?{ 3 M^ + 2 M‘^{yipi2 + ITpis 
+ 17p23 - 5Qi - 5Q2 - 4Q3) + 2(18pi2 + 13pi3 + 13p23 - 7p23Ql -Ql- 8 p 23 
XQ2 — 2Q1Q2 — Q2+ Pl2(31pi3 + 31p23 ~ llQl ~ IIQ2 ~ I2Q3) — JP23Q3 
— I4Q1Q3 — 14 :Q 2 Q 3 — 2Q3 + Pl3(18p23 ~ 8Q1 — 7 {Q 2 + Q3)))) 
+ 3 ieaf 3 ^SPtP 2 P 3 Q^iPl 3 “ P 23 - Ql + Q2)) , (5.85) 

2048mM(3m® + 2 pl^p 23 + 4pi3P23 + 2pL + 5pi2P23Qi + 7 pi 3 P 23 Qi + QplsQi 
-6P23Q? - QPI2Q2 - ‘2IP12P13Q2 - I5P13Q2 - 8P12P23Q2 - I8P13P23Q2 - 4P23 
xQa + 6P12Q1Q2 + J‘ 2 pi 3 QiQ 2 ~ 7P23Q1Q2 + 6P12Q2 + I5P13Q2 + 2P23Q2 
+ M^(3p23Ql + Pl 2 {p 23 ~ 8Q2) — Pl 3 {p 23 + 6Q2 “ 3Q3)) + 3p^2Q3 + QpuPnQs 
+6P13Q3 + 6P12P23Q3 + 4P13P23Q3 ~ 6P13Q1Q3 — 5P23Q1Q3 + 3P12Q2Q3 + 9 pi 3 
XQ2Q3 ~ 8P12Q3 — 6P13Q3 + + 2pi2 + 6pi3 + 9p23 + llQi ~ 23Q2 

+ I5Q3) + m^(pi2 + 3pi3 + 10pi3P23 + 8P23 + + l 7 p 23 Ql - QQl - 38pi3 

XQ2 — 2IP23Q2 + 7 )QiQ2 + I5Q2 + 77P13Q3 + 9^23^3 ~ I7Q1Q3 + 2Q2Q3 

— 11 Q\ + M^(pi 2 — pi 3 — P23 + 3 Qi — 6Q2 + 6(53) + Pi 2 ( 2 pi 3 + 2 p 23 + 

-23Q2 + I6Q3)) - 3 iea^^ 5 p‘lp 2 PlQ\ 2 m‘^ + pia + P13 + 2p23 - Q2 - Qs)), 

(5.86) 

6144(3m^ - P12 - 2pi2(pi3 + P23) + (pia + P23Y + 4 m^{-pi 2 + P13 + P23 )), 

(5.87) 
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fSV,2332 _ ^5\/,1441| ^5F,3443 _ fSV,122l\ 

J LLRR ~ J LLRR \pi-^V2 ’ J LLRR ~ J LLRR lpi-H-p2,mi->-m • 


( 5 . 88 ) 
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